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The  view,  opinions,  and/or  findings  contained  m 
this  report  are  those  of  the  authors  and  should 
not  be  contrued  as  an  official  Department  of  the 
Army  position,  policy,  or  desision,  unless  so 
designated  by  other  documentation. 


Abstract  and  Summary 

This  report  describes  and  evaluates  a  new  way;  to 
combine  the  two  interference-resistant  techniques  of  pseudo¬ 
noise  (PN)  spread  spectrum  and  adaptive  (antenna)  arrays.  Each 
element  in  the  array  feeds  a  separate  PN  matched  filter,  and 
the  sampled-matrix-inversion  (SMI)  algorithm  is  used  on  the 
matched-filter  outputs,  along  with  multiple  steering  vectors 
corresponding  to  apriori  beams.  The  SMI  calculation  is  done 
on  a  sliding-window  basis,  where  the  noise  covariance  matrix 
is  estimated  using  matched  filter  output  samples  that  precede 
the  "current-sample,"  and  the  resulting  multiple  optimum 
weights,  computed  by  matrix  inversion  multiplied  by  the  multiple 
steering  vectors,  are  applied  to  the  current  sample.  The 
number  of  samples  used  in  the  covariance  matrix  estimation  is 
at  least  twice  the  number  of  antenna  elements. 

The  evaluation  of  this  PN  matched  filter  plus  adaptive 
array  technique  was  done  by  using  a  4-element  linear  array  and:  j 
1)  straightforward  matrix  calculations;  2)  a  computer  computa-  j 
tion  of  output  signal- to-noise  ratio;  SNRq,  versus  angle  when 
the  theoretically  :ect  noise  covariance  matrix  was  used;  and 
3)  a  computer  simulc.v  m/demonstration  that  estimated  the  noise 
covariance  matrix  from  computer-generated  interference  samples.  , 
It  was  found  that,  as  expected,  the  separation  of  j 


interference-only  samples  from  interference-plus-signal  samples 


afforded  by  the  pulse-compression  nature  of  the  matched  filter, 
permits  full  combined  processing  gain  (PN  plus  adaptive  array) 
from  the  outset,  in  the  acquisition  mode.  We  confirmed  that 
the  noise  covariance  matrix  estimate  requires  using  only  sample 
sizes  on  the  order  of  twice  the  number  of  elements  for  3  db  loss 
from  optimum. 

We  used,  as  a  theoretical  performance  measure,  the 
output  SNRq  as  a  function  of  angle,  assuming  the  steering  vector 
points  in  the  direction  of  the  "current  angle."  Realizing  this 
performance  in  practice  would  require  a  (slightly)  different 
steering  vector  for  each  (say)  degree.  In  the  effort  here  we 
demonstrate  that  one  need  only  use  about  4  steering  vectors  to 
achieve  near-optimum  SNRq  for  the  entire  range  of  angles. 

Thus,  a  near-optimum  combined  processing  has  been 
effected  which  can  be  viewed  as  sliding-window  in  both  the 
time  axis  and  in  the  spatial  domain  (due  to  the  multiple  steering 
vectors).  The  theoretical  computer  computations,  which  used  the 
ideal  covariance  matrix,  were  used  to  portray  SNRQ  versus  angle 
for  a  variety  of  single  and  multiple  interferer  conditions,  for 
a  single  signal.  The  computer  simulation/demonstration  docu¬ 
mented  the  performance  for  two  single-interferer  examples,  and 
one  two-interferer  example. 

These  results  indicate  that  using  the  method  evolved 
here,  along  with  the  4  parallel  steering  vectors,  may  likely  be 
the  best  way  to  combine  PN  waveform  processing  and  adaptive  array 
processing.  The  advantages  are:  1)  the  combination  does  not 


initially  (before  sync)  attempt  to  null  the  signal,  but  offers 
immediately  the  full  combined  processing  gain;  2)  the  response 
to  either  signal  turn-on  or  interferer  variation  is  near  instan¬ 
taneous,  and  does  not  depend  on  loop  time-constants;  3)  there  is 
no  constraining  relation  between  interferer  size  and  time-response, 
as  there  is  for  a  closed-loop  system,  and  4)  there  is  no  "great 
race"  problem  between  modem  synchronization  and  array  nulling, 
as  there  is  when  a  PN  reference  is  used  in  a  closed-loop;  and 
5)  the  technique  here  should  be  more  resistant  to  specific  and 
determined  interferer  attacks  (such  as  blinking) .  It  is  recom¬ 
mended  that  this  method  be  further  pursued,  to  further  develop 
and  document  its  advantages. 
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Fig.  49.  SNR  versus  angle,  when  covariance  matrix  is  estimated  from  arriving  signals, 
with  five  Interferers  of  Nj  *  10^W  at  70°,  120°,  190°,  230°,  and  290  ,  for 
averaging  over  2,4,8,16,32,  and  200  samples.  (Na  a  0.1W,  S  *  1W  at  1  ) 
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PN  MATCHED  FILTERS  IN  ADAPTIVE  ARRAYS 


1.  Introduction 

The  overall  objective  in  this  effort  is  to  achieve 
more  effective  and  robust  arrangements  for  combining  spread 
spectrum  and  spatial  processing  gain  in  ECCM  communication 
systems.  We  assume  here  that  the  spread  spectrum  is  pseudo¬ 
noise  (PN) ,  and  investigate  the  potential  of  using  identical 
PN  matched  filters  at  each  antenna  element,  followed  by  a 
spatial  nulling  algorithm.  Thus,  a  significant  feature  of 
the  approach  here  is  that  the  waveform  processing,  and  hence, 
temporal  processing  gain,  precede  the  spatial  processing. 

The  spatial  processing  most  suitable  for  this  approach  is  some 
form  of  direct  matrix  inversion?  we  have  used  the  straight¬ 
forward  sampled  matrix  inversion  (SMI)  for  this  effort. 

2.  Background 

Two  techniques  for  providing  interference'*" -resistant 
communications  against  either  unintentional  or  intentional 
interference  are  spectrum  spreading  and  antenna  gain.  The 
antenna  gain  here  refers  to  the  gain  in  the  desired-signal 
direction  viz-a-viz  the  gain  in  the  interferer  direction.  In 
many  communication  applications  either  the  terminals  themselves 

1.  Interference  refers  to  any  man-made  electro-magnetic  signal, 
as  distinguished  from  nature's  ubiquitous  Gaussian  noise 
(WGN) .  It  is  any  signal  from  a  point  source  (plane  wave) 
which  is  not  correlated  with  the  desired  PN  signal. 


2. 


or  the  interferers  move  in  physical  space  in  an  unpredictable 
manner,  so  that  the  antenna  gain  must  be  adaptive,  using 
adaptive  arrays. 

The  essential  idea  in  spectrum- spreading  is  to  dis¬ 
perse  the  intended  signal  over  a  frequency-space  that  is 
large  relative  to  that  required  by  the  information  data  rate, 
which  permits  a  temporal  signal  processing  advantage  against 
any  interference  (jamming  or  RFI)  that  is  subject  to  elec¬ 
tronic  amplification  limitations^-  and  is  not  correlated  with 
the  PN  signal.  The  adaptive  array  reduces  the  antenna  gain 
in  the  direction  of  the  interference  in  such  a  way  as  to  maxi¬ 
mize  the  SNR  for  the  desired  signal. 

Past  interference-resistant  communications  have 
usually  used  spread-spectrum  (utilizing  waveform  or  temporal 
processing) .  More  recently,  adaptive  antenna  techniques  (utili¬ 
zing  spatial  processing)  have  been  combined  with  spread  spectrum 
to  provide  increased,  and  more  robust,  interference-resistance. 
Combining  the  two  techniques  offers  the  potential  of  realizing 
the  sum  (in  db)  of  the  two  processing  gains  against  either 
the  intentional  or  unintentional  interferences. 

If  frequency-hopping  (FH)  is  the  spread-spectrum 
technique,  then  the  antenna  array  is  "adapted"  in  a  step-by- 
step  fashion,  slightly  ahead  of  the  frequency-hop  path.  If 
pseudo-noise  (PN)  is  the  spread  spectrum  technique,  then  the 

1.  There  is  no  temporal  processing  gain  against  white  Gaussion 
noise  (WGN) ,  such  as  receiver  noise,  because  the  noise  is 
already  present  at  all  frequencies.  Increasing  signal  band¬ 
width  then  also  increases  the  total  noise  received  in  the 
signal  bandwidth.  Any  interference  subject  to  electronic 
amplication  is  limited  by  the  power  dissipation  capability 
of  electronic  devices;  in  this  case  the  total  interference 
from  any  given  source  is  fixed.  Increasing  signal  bandwidth 
then  allows  the  intended  receiver  to  reject  a  longer  fraction 
of  the  interference  power  from  such  a  source. 


3. 


antenna  is  adapted  on  either  a  time-continuous  or  a  sampled 
basis  corresponding  to  the  symbol  rate  of  the  digital  trans¬ 
mission. 

The  fundamental  operations  of  an  adaptive  antenna 
(alone)  are  depicted  in  Fig.  1.  Signals  received  in  each 


element  channel  are  weighted  by  a  complex  coefficient  and 
summed  together  to  form  the  output  of  the  array.  Temporarily 
assume  that  only  interference1  sources  are  present,  so  that  any 
ensuing  receiver  deals  with  directional  interference  power  (which 
is  to  be  nulled)  plus  ordinary  receiver  thermal  noise.  It 
is  useful  to  describe  the  desired  element  weighting  (complex) 
vector  in  terms  of  the  "covariance  matrix"  of  the  impinging 
noise  and  interfering  signals,  plus  a  possible  steering-vector. 

1.  Although  many  adaptive  arrays  are  intended  to  operate  against 
jammers,  the  algorithms  adapted  work  equally-well  against  un¬ 
intentional  interference  that  may  be  present,  including  other 
friendly  users.  Therefore,  we  will  use  the  general  term  inter 
ference,  rather  than  jammer,  throughout  this  report. 


4. 


It  is  the  intent  of  the  covariance  matrix  to  des¬ 
cribe  in  systematic  form  all  of  the  information  about  the 
interference  environment  (see  Ref.  1.)  Let  the  covariance 
matrix  M  have  elements  (m^) 

M  =  (nijtl)  (1) 

The  element  value,  m^,  is  the  expected  value  of  the 
product  of  the  (complex  envelope  representation)  noise  re¬ 
ceived  (jamming  plus  receiver  noise)  by  the  kth  and  1th 
antenna  elements: 

* 

mj^  =  E  (mjcm1)  (2) 

This  implements  the  definition  of  covariance  (see 
any  communications  textbook) ,  and  physically  represents  the 
crosspower  in  the  envelopes  of  the  two  interference  signals. 
Thus,  m^  is  the  cross-power  in  the  envelopes  of  antenna 
elements  k  and  element  1.  These  element  values  clearly  con¬ 
tain  directional  interference  information  since  any  RF  delay- 
pattern  set  up  by  the  point-source  and  antenna  array  con¬ 
figuration  determine  the  averaged  cross-product  of  each 
element-pair.  In  these  terms,  a  weight-vector  which  will 
theoretically  yield  the  best  signal-to-interference  ratio 
for  a  single  (desired)  signal  case  is  given  by: 

MW  =11  S*  (3) 

where:  M  =  covariance  matrix  of  the  interference 

Vi  *  arbitrary  nonzero  complex  constant 

* 

S  =  complex  conjugate  of  the  desired-signal  vector 
(from  the  array),  i.e.,  a  steering  vector. 

W  «  complex  weight  vector 


For  radar  cases,  the  S  vector  for  electronically- 
steered  arrays  is  essentially  a  beam-steering-vector .  For 
communication  cases,  where  the  direction  of  the  desired 
signal  is  not  known  a  priori,  the  S  vector  may  either  cor¬ 
respond  to  an  omnidirectional  pattern,  or  to  a  hypothesized 


desired  signal  direction.  In  either  case,  the  optimum  weight 
vector  can  be  written  asA 

W  =  y  M_1S*  (4) 

It  is  useful  to  describe  any  adaptive  antenna 
processing  in  terms  of  its  handling  of  the  matrix-inversion 
contained  in  Eq.  4.  It  is  fundamental  to  distinguish  between 
open-loop  and  closed- loop  implementation  approaches,  which 
are  depicted  in  Fig.  2.  Two  basic  approaches  are:  1)  imple¬ 
ment  Eq.  4  directly  by  sampling  the  observed  M  matrix,  and 
inverting  the  observed  matrix  (shown  in  Fig.  2(a));  2)  form 
a  servomechanism-like  closed  loop,  and  use  a  processing- 
algorithm  in  that  loop  (shown  in  Fig.  2(b)). 

The  open-loop  technique  is  called  sampled-matrix 
inversion  (SMI),  and  is,  of  course,  a  digital  processing 
technique.  This  method,  since  it  is  open-loop,  requires 
relatively  high  arithmetic  precision  to  obtain  good  results. 

To  obtain  an  accurate  estimate  of  the  covariance  matrix, 
it  has  been  found  (Ref.  9)  that  at  least  2K  such  sample  ma¬ 
trices  should  be  calculated  and  averaged  for  the  estimated 
covariance  matrix.  (K  is  the  number  of  elements.)  Also,  the 
covariance  samples  should  be  obtained  when  no  desired  signals 

T.  M"1  is  the  inverse  of  M.  It  will  exist  if  ambient  noise  is  j 
present.  It  is  often  close  to  singular  with  strong  interferend 
2.  This  assumes  a  known  steering  vector.  In  actual  operation  a  ] 
somewhat  larger  number  of  samples  may  be  needed  (Ref.  11) .  j 
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Fig.  2.  Adaptive  Receiving  Antennas:  a)  open-loop; 


are  present.  Otherwise  the  array  gain  in  the  desired  signal 
directions  will  also  be  reduced. 


The  advantages  of  the  algorithm  are:  (1)  It  gives 
a  clearly  specified  algorithm  which  is  optimal  for  any  assumed 
noise  environment;  (2)  the  response  to  a  change  in  noise  en¬ 
vironment  can  be  very  fast.  The  fundamental  limitation  in 
response  time  is  the  time  required  to  obtain  the  independent 
samples  for  the  covariance  matrix. 

The  closed-loop  (servo-like)  arrangement  in  Fig.  2(b) 
has  been  popular  in  radar  work,  and  has  been  the  starting 
point  for  communication  applications.  Such  a  closed-loop 
arrangement  has  the  self-adaptive  properties  inherent  in  any 
servo-like  loop.  It  can  be  shown  (Ref.  1)  that  the  feed¬ 
back  rule  that  maximizes  the  signal-to-interference  ratio  is 
the  "least  mean  square  (LMS ) "  algorithm.  In  effect,  the 
loop  attempts  to  adapt  to  the  situation  where  the  feedback 
"error"  signal  is  zero.  Physically,  the  loop  attempts  to 
place  a  (relative)  null  on  every  point  source  in  its  field 
of  view  by  adjusting  the  antenna  weights  so  that  the  signal 
output,  after  combining,  is  a  least-mean-square  best  fit  to  a 
stored  copy  of  the  expected  signal  waveform  (zero,  if  ambient 
noise  is  neglected).  The  stronger  the  signal,  the  deeper  the 
null.  Such  a  (loop)  antenna  continuously  adjusts  its  own 
pattern  by  means  of  feedback  as  it  operates. 

It  is  interesting  to  note  that  the  asymptotic  values 
of  the  complex  weights,  with  the  LMS  algorithm,  are  exactly 
the  covariance  matrix  (M)  values  given  in  Eq.  4.  Thus  the 
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theoretical  underpinnings  of  the  open-loop  (SMI)  and  the 
closed- loop  (LMS)  are  the  same.  However,  their  actual 
implementation,  especially  when  combining  with  waveform 
(spread  spectrum)  processing,  and  their  response  to  jammers 
differ  substantially  (as  we  shall  see)  .  Also  we  may  note 
that  both  implementations  attempt  to  "spatially  whiten"; 
that  is,  reduce  gain  the  direction  of  point  sources  so  that 
all  directions  appear  to  have  a  noise  level  equal  to  the 
ambient  noise  level. 

Thus  far,  the  closed-loop  version  has  been  pursued 
when  combining  PN  spread  spectrum  with  adaptive  antennas. 

A  procedure  of  particular  interest  has  been  using  the  "Widrow 
LMS  processor"  with  a  PN  signal  as  the  reference  signal.  Fig.  3 
shows  a  basic  arrangement  (Ref.  2)  where  a  PN  signal  forms  the 
reference  signal.  When  the  PN  code  is  synchronized,  the  de¬ 
sired  signal  part  of  the  array  output  passes  through  the  code¬ 
reference  loop  without  change,  but  the  modem  data  can  be  ex¬ 
tracted  as  shown.  The  first  mixer  removes  the  PN  code  so  that 
the  signal  will  pass  through  the  data  bandwidth  filter,  the 
limiter  fixes  the  amplitude,  and  the  second  mixer  puts  the 
code  back  on.  In  this  way,  the  PN  code  prevents  the  adaptive 
array  from  attempting  to  null  the  intended  signal.  Rather, 
nulls  are  placed  on  any  interference  source  which  are  uncor¬ 
related  with  the  reference  signal.  Note  that  the  PN  modem 
involved  is  of  the  active-correlator  type.  In  Ref.  3  the  same 
basic  technique  is  used  except  that  now  multiple  delays  are 
used  with  each  antenna  element  so  as  to  implement  better  wide¬ 
band  null  steering. 


Using  LMS  Processor 


Two  significant  aspects  of  the  arrangement  in  Fig.  3 


are:  1)  the  spatial  processing  precedes  the  waveform  of 

temporal  (spread  spectrum)  processing;  and  2)  the  architecture 
is  that  of  a  closed  loop,  with  the  PN  modem  arranged  to  form 
the  "known  signal"  reference  for  the  loop  so  that  the  intended 
signal  is  not  nulled. 

An  alternative  to  the  closed-loop  approach  of  Fig.  3 
was  suggested  in  Ref.  4,  and  uses  an  active  PN  correlator  ahead 
of  the  adaptive  array  processing,  as  depicted  in  Fig.  4.  This 
arrangement,  in  its  basic  form,  is  an  open-loop  technique  and 
hence  suggests  the  use  of  direct  Sampled  Matrix  Inversion  (SMI) 
As  noted,  significant  aspects  of  any  such  open-loop  technique 
are:  1)  it  is  inherently  digital,  requiring  both  sampling  and 
quantizing;  2)  it  can  react  instantaneously  as  opposed  to 
exhibiting  a  transient  response. 


Fig.  4.  Open  loop  PN  technique  using  Direct  Samples  Matrix 
Inversion 
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A  recent  report'*' *  (Ref.  5)  includes  a  steady-state 
investigation  of  the  arrangement  shown  in  Fig.  4,  and  uses 
the  descriptive  terms  "code-stripping"  to  describe  this 
arrangement.  We  will  discuss  this  further  in  the  next  section. 

In  Ref.  6  the  SMI  algorithm  is  evaluated  for  use 
with  a  frequency-hopping  spread  spectrum  temporal  processing. 

In  that  case  the  receiver  can  estimate  the  noise  covariance 
matrix  one  step  ahead  of  the  frequency-hopping  pattern,  and 
apply  the  calculated  weights  of  the  signal-present  subse¬ 
quent  interval.  As  noted,  the  interest  here  is  in  cases 
where  PN  is  the  spread  spectrum  technique. 

3.  PN  Matched  Filters  Preceding  SMI  Array  Processing 

In  Ref.  7  we  proposed  using  PN  matched  filters  at 
each  antenna  element,  and  envisaged  using  multiple  (time- 
delayed)  samples  at  the  output  of  each  matched  filter  to 
feed  a  (then  as  yet)  undetermined  null-forming  and  beam¬ 
forming  processing.  In  the  early  stages  of  this  effort  we 
realized  that  direct  SMI  was  an  appropriate  spatial  process¬ 
ing  for  processing  samples  of  the  matched  filter  outputs. 

A  filter  matched  to  a  maximal  length  (see  Ref.  8) 

PN  sequence  has  the  unique  feature  that,  assuming  that  two 

2 

full  periods  cf  ~he  sequence  are  transmitted,  the  output  ex¬ 
hibits  ideal  "pulse  compression."  That  is,  all  the  energy 
from  a  single  period  of  the  received  (noiseless)  sequence 
appears  at  a  single  peak,  and  the  preceding  samples  values 

1.  The  work  reported  in  Ref.  5  was  conducted  during  the  same 
period  as  the  work  of  this  report. 

One  can  avoid  the  need  for  2  periods,  using  specially- 
constructed  matched  filters.  This  implementation  aspect 
is  beyond  the  scope  of  this  study. 
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AMPLITUDE  OP  SAMPLES 


(at  independent  times)  can  be  arranged  to  be  zero.  The 
idealized  time  description  of  sampled  values  of  the 
matched  filter  output,  assuming  no  noise  is  present,  is 
sketched  in  Fig.  5.  Note  that  if  the  phase  modulation  used 
is  slightly  less  than  180°,  the  signal  output  will  be  exactly 
zero  between  signal  peaks. 
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Fig.  5.  Sampled  Envelope  Values  of  Signal  Autocorrelation 


This  pulse  compression  feature  nf  the  matched 
filter  suggests  that  the  SMI  algorithm  is  sensibly  applied  to 
the  K  matched  filter  outputs  in  the  following  way.  Assume 
momentarily  that  the  (signal)  matched  position  is  the  "current 
sample"  value.  Then  the  immediately  preceding  (to  the  left) 
independent  samples  must  contain  noise-only.  About  2K  preceding 
samples  from  each  of  the  K  matched  filter  outputs  are  used  to 
estimate  the  noise  covariance  matrix.  The  optimum  weights  are 
computed  by  inverting  this  matrix  and  multiplying  by  an  apriori 
beam  steering  vector,  to  be  described.  These  K  complex  values 
are  then  used  to  weight  the  K  "current-sample"  values.  The 


process  just  described  can  be  operated  in  a  sliding-window 
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fashion,  so  that  this  procedure  is  suitable  for  the  acquisition 
mode,  where  the  filter  "matched  position"  is  not  known  apriori 
but  must  be  determined  by  an  ensuing  threshold-crossing  algorithm. 

The  significant  point  to  be  made  here  is  that  the 
pulse  compression  nature  of  the  matched  filter  output  permits 
separating  signal  plus  noise  samples  from  noise-only  samples, 
while  still  allowing  use  of  a  continuous  signal.  Furthermore, 
the  noise  samples  used  are  a  subset  of  those  that  occured 
during  the  signal-epoch  corresponding  to  the  matched  filter 
length. 

Thus  the  algorithm  pursued  in  this  effort  is  pictured 
schematically  in  Fig.  6.  In  addition  to  using  a  sliding- 
window  in  the  time  axis  (or  temporal  domain) ,  we  use  an 
equivalent  feature  in  the  spatial  domain  consisting  of  using 
multiple  apriori  beams  (or  steering  vectors).  These  apriori 
beams  mean  that,  in  the  absence  of  any  point-source  interferers, 
multiple  beams  will  be  formed  with  the  multiple  elements  of  the 
array.  The  number  and  orientation  of  such  beams  are  decided  in 
advance,  and  do  not  vary  with  the  interference  scenario.  When 
combined  with  adaptive  processing,  the  use  of  multiple  apriori 
beams  amounts  to  performing  multiple  weight-calculations,  which 
results  in  multiple  output  time-tracks.  That  patticular  process¬ 
ing  which  has  a  beam  placed  in  the  general  direction  of : the  in¬ 
tended  signal  should  experience  the  highest  spatial  processing 
gain,  from  among  the  multiple  processsings. 

Thus,  with  the  understanding  that  the  samples  used 
for  the  covariance  matrix  are  separate  from  the  samples  to 
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which  the  computed  weights  are  applied,  we  can  write: 


where :  W 
-!■ 

M 


S 

k  SL 


Kt  -  (A>'1  skl  (5) 

the  estimated  optimum  weight  vector  for  the 
Jlth  apriori-beam 

K  by  K  covariance  matrix,  estimated  by  2K 

samples  from  +  .  to  t  . 

n-l  n-2K-i 

a  K  by  SL  matrix  that  represents  £  different 
apriori  beams  or  "beamformers. " 


If  one  multiplies  the  current  samples,  t^,  by  the 
l  weight  vectors,  one  obtains  SL  outputs,  ,  which  is  a  vector 
having  one  element  for  each  apriori  beam: 


where : 


h.  ~ 

£th  output  of  SMI  processing 
output  of  K  matched  filters 


(6) 


4.  Theoretical  Performance  Evaluation 

Here  we  assume  a  linear  array  of  K  elements,  with 
element  spacings  d,  as  sketched  in  Fig.  7. 


Fig.  7.  Linear  Array  of  K  Elements. 
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4 . 1  Array  Performance  Measure 

To  begin,  we  consider  only  one  interferer,  at 
electrical  angle  0  .  The  array  performance  measure  used  here 
is  the  output  signal  to  noise  ratio  when  the  beamformer  is 
pointed  in  the  direction  of  the  unit  signal,  at  angle  ©  . 

We  normalize  so  that  the  total  received  ambient  noise  is 

1. 

With  the  single  interferer,  a  normalized  theoreti¬ 
cal  covariance  matrix  can  be  shown  (see  App.  A)  to  be: 


M  =  NaI  +  NiRq 


(7) 


where : 


N  = 


I 

N, 


ambient  (mean-square)  power  per  antenna  element, 
identity  matrix 

interference  (mean  square)  power  per  antenna  element. 
*T 


vnv*T  -  J_  A  (Rt't'-co'-'ie 
0  0  K 

1  I  1  I 


i 


’0 


1 


0  = 


2irfd  ! 


-3© 


;-zj0 


-(K-i)  jel 


electrical  phase  shift  in  radians  between 


Si 


c  cosa 

elements  separated  by  a  distance  d,  for  a 
ane  wave  arriving  at  an  angle  a'  degrees  with  respect  to  the  line 
the  array. 


Then  the  inverted  covariance  matrix  is: 

.-1  _ 


M 


».(I  '  V"i  R®) 


(7a) 


C 


It  is  germane  to  ask  what  interferer  bandwidth 
implications  arise  from  writing  the  receiver  voltage  vector, 
VQ  ,  in  terms  of  a  plane-wave  at  a  given  RF  frequency.  This 
technique  essentially  assumes  that  the  envelope  of  the 
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arriving  interference  is  the  same,  for  all  elements  of  the 
array.  This  will  certainly  be  true  for  any  narrowband  inter- 
ferer,  and  will  be  approximately  true  for  random  PN-coded 
interferers  and  Gaussian  matched-spectrum  interferers.  In 
the  following,  we  will  be  implying  that  the  interferers  are 
either  random-coded  PN  or  Gaussian  matched- spectrum. 

Note  that  a  random-coded  PN  interference  signal 
can  be  treated,  for  many  purposes,  as  approximately  a  piecewise 
CW  signal,  since,  for  the  vast  majority  of  the  time  axis,  each 
element  in  the  array  is  seeing  the  same  "chip." 

From  classical  theory  (see  App.  A) ,  and  using  the 
above-described  procedure  for  forming  multiple  weights  by 
using  multiple  steering  vectors,  the  weights  for  a  steering 
angle,  o',  are: 

V  =  (M~1)  V  (8) 

where:  W  »  =  optimum  weights t using  an  apriori  beam  in 

0  the  direction  0  .  f 

Vq'  =  same  vector  as  before,  except  0  replaces  ® 

The  resulting  normalized  beam  point  pattern  for  a  signal  from 
direction  <p  is: 

*T 

»e  %  <9> 

where:  <p  =  2  cos  {5 

g  =  any  physical  angle,  in  degrees,  between  the  signal 
source  and  the  line  of  the  array. 

Since  we  have  chosen  to  evaluate  the  performance,  in  any  given 
direction  g ,  as  the  magnitude  of  the  beam-pattern  when  the 
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apriori  beam  (steering  vector)  points  in  that  direction 
(the  0*  =  <j>  )  ,  the  output,  power  due  to  the  signal  in  that  same 

direction  is: 

*T  *T  *T  *T 

Po  <j>  =  M'\  V,  =  V*  M"1  V*  =  V*  M'\  (10) 

where:  PQ(j)  =  output  power,  due  to  interference,  in  direction  <t> . 

We  can  now  note  that  the  output  SNR  equals  P  ^ since:  1)  the 

2 

output  signal  power,  for  electrical  direction  <j> ,  is  P  .  , 
and  2)  the  output  noise  power  for  electrical  direction  i)>  is 
Po^.  Hence,  the  on-beam  output  SNR  is  simply  P  ^  . 


Equation  (10)  implies  that  the  noise  output  power, 
for  a  given  angle,  is  not  constant,  but  is  a  function  of  the 
steering  direction.  To  achieve  the  same  SNR  gain,  but  have 


spatial  whitening  (which  is  a  superior  model  concept)  one  lets: 


-1  * 
(M  XV  ) 


W  = 


-1  * 
(M  V  ) 


*T  -1 
V*  XM  % 


Then,  let  Z  represent  the  output: 


Z  «  W^Ty  (12) 

In  noise  alone: 

E  {  j  Zj  2  }  *  1  (13) 

for  all  straight  angles, 
and  in  unit  signal,  on-beam: 


1 


(M  xVd>)  V* 

E{  Z  }  = 

yi»  p 

% 

(14) 

and 

|E(Z)  |2  = 

Po$ 

(15) 
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This  is  the  output  signal  for  unit  input  power  per  element. 
Since  output  noise  is  normalized  to  1,  it  is  also  the  output 
SNR  when  the  steering  vector  is  in  the  direction  of  the  desired 


signal. 


In  the  case  of  a  single  interference,  in  addition  to 


ambient,  the  P  ,  can  be  reduced  to  a  single  equation  (see 


App.  A),  and  appears  as: 


.  i_  _*i _ r 

N  "  N  (N+nT  ' 
a  a  a  1  J 


sin  Kir- —  (cos  a  -  cosa 

” — : — hr~. - 

K  sm  it  —  (cos  a  -  cosa 
c 


ii 

) 

T 


For  a  signal  power,  S,  per  element  and  an  array  of  K  elements, 


this  becomes: 


=  K 


/ s  \  r,  Si/S  f sln  Kt  5~ 

\y  V  “a  U 

*—  f  r* 


—  (cosa  -cosa  ) 

C  T 


(cosa  -cosa^.  ) 


(17a) 


where:  rr  =  Received  ambient  —  per  element  (omnidirectional) 

Na  Na 

N 

I  =  Received  directional  interference  to  signal  ratio 
c 

per  antenna  element  (omnidirectional) 

a  =  angle  between  line  to  interference  and  line  of  the 
array1  in  radians 

a  =  angle  between  line  to  desired  signal  and  line  of 

x 

the  array  in  radians 


The  squared  term  on  the  right  will  be  recognized  as  the  sine 
periodic  sampling  function  which  is  small  for  all  angles  except 
when  the  denominator  is  zero.  Then  sine  *=  1  and  the  output 

1.  Note:  This  differs  from  the  more  usual  case  where  angles 
are  measured  from  the  array  normal,  a  *  “  7  where  Y  is 

the  angle  from  normal. 
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C 

becomes  SNRoufc  =  K  -■  - ~N  ■■  .  This  happens,  for  example,  whenever 

I  a 

cosa  =  cosa^  .  One  example  is  when  the  signal  and  interference 
come  from  the  same  direction  (a  =  a^).  The  interference  is  sup¬ 
pressed  completely  whenever  sin  (Kit  cos  a  -  cos  a^)  =0,  while 
sin  ir  —  (cosa  -  cos  a  )  0.  The  important  feature  of  the 

equation  is  that  in  practical  cases,  the  interference  is  almost 
completely  suppressed  almost  everywhere  except  in  the  immediate 
direction  of  the  interference. 

A  series  of  selected  examples  for  this  single-interferer 
case  are  collected  in  Ref.  12. 

In  general,  one  can  note  that: 

M  =  E(YY*)  =  ^  |in(t)|2  V  (0n)  V*  (0n)  +  NaIR  (17) 

n 

where  we  have  assumed  independence  of  interferers  so  that  all 
expected  values  for  cross-product  terms  are  zero. 

Where:  i  (t)  =  interference  time  waveform  from  n^*1  interferer 

n 

V  (0  )  =  voltage  vector  received,  at  angle  0  ,  from 

n  th  .  .  _  n 

n  interferer 

N  Iv  =  covariance  matrix  of  ambient  noise, 
d 

We  found  no  single  expression  which  can  relate  theoreti¬ 
cal  SNR  performance  to  angle  for  the  multiple  interferer  case, 
so  we  resorted  to  using  a  computer  computation,  the  results  of 
which  is  now  described. 

4.2  Theoretical  Computer  Computations 

A  computer  program,  in  BASIC  language,  was  written 
to  be  run  on  our  own  CEL  Line  8  minicomputer.  The  intent  of  this 
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program  was  to  find  the  theoretical  predicted  results  for  the 
algorithm  and  linear  arr  being  used  here  for  general  cases, 
where  more  than  one  interferer  is  present. 

Therefore,  the  theoretical  covariance  matrix  was 
used;  this  was  formed,  after  the  interferer  directions  and  powers 
are  stated,  by  deterministrative  computation,  using  the  usual 
trigonometric  relations. 

The  first  program,  CMAT4,  assumes  that  an  infinite  number  of 
steering  vectors  are  available,  since  the  SNR  versus  angle  4> 
always  assumes  that  the  steering  vector  is  in  the  direction 
of  the  intended  signal  no  matter  what  signal  direction  is  chosen. 
A  four-element  half-wavelength  spaced  linear  array  is  assumed. 

The  user  specifies  the  power  and  direction  of  each  directional 
interference,  and  the  ambient  power.  The  matrix  inversion  is 
accomplished  by  using  elementary  row  operations.  The  CMAT4 
computes  the  dot-product  performance  measure  for  0  to  180 
degrees  in  10  degree  steps,  then  allows  the  user  the  option  of 
selecting  a  different  signal  direction  range  or  stop,  or  re¬ 
cycling  for  a  new  interference  description. 

It  must  be  emphasized  that,  in  the  results  here,  the 
covariance  matrix  was  not  obtained  by  an  estimate  using  the 
incoming  samples,  but  rather  was  formed  by  initially  stating  the 
interferer  directions  and  powers,  and  then  using  trigometric 
plane-wave  relations  to  deterministically  find  the  covariance 
matrix. 
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Note  that  an  operational  communication  system 
does  not  know  either  the  noise covariance  matrix  or  the 
exact  direction  from  which  a  desired  signal  is,  or  will 
be  arriving.  It  must  estimate  the  covariance  matrix  from  a 
set  of  noise  samples.  It  can  choose  that  one  of  the  set  of 
steering  vectors  which  maximizes  the  SNR  for  its  approximate  co- 
variance  matrix,  but  this  generally  is  not  exactly  the  correct 
steering  vector  for  the  actual  signal. 

In  this  section,  we  compute  the  exact  optimum  SNR 
based  on  the  known  parameters.  In  the  next  sections,  we  will 
estimate  how  close  an  actual  system  working  with  finite 
numbers  of  noise  samples  can  come  to  this  optimum  performance. 

Here  we  discuss  the  SNR  performance  data  for  numerous 
examples.  In  the  next  section  we  discuss  the  antennae  power 
patterns  which  are  obtainable  from  the  programs.  The  SNR 
in  db  is  displayed  as  a  function  of  signal  direction  for 
the  optimum  case  where  the  steering  angle  is  pointed  in  the 
signal  direction.  As  noted  in  Sec.  4.1,  the  analysis  and  com¬ 
putation  here  applies  accurately  to  narrowband  interferer 
cases,  and  approximately  to  both  random-coded  PN  and  Gaussian 
spectrum-matched  interferers. 

On  these  curves  S/N_  is  the  S/N  observed  at  the 

®  cl 

output  of  each  of  the  4  matched  PN  filters  when  only  ambient 
noise  is  present.  NJ/S  is  the  NI/g  observed  at  the  output  of 
each  of  the  matched  PN  filters.  The  output  (S/N)0  *s 
the  actual  output  signal  to  noise  ratio. 
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In  Fig.  8  the  ambient  noise  at  the  output  of  each 

s 

matched  PN  filter  is  assumed  to  give  an  ambient  —  of  10  db. 

Na 

In  the  absence  of  all  jammers  the  performance  would  be  +16  db 
in  all  receiving  directions,  as  shown  by  the  dotted  line.  The 
presence  of  a  30  db  Nj/S  interferer  at  +45°  (measured  counter¬ 
clockwise  from  the  +  x  axis),  results  in  the  decrease  in 
shown  for  signals  near  +45°.  Since  the  ambient  is  much 
smaller  than  the  interferer,  the  minium  (S/N)q  is  essentially 
that  due  to  the  interferer,  -30  db.  Example:  S  =  1  Watt 


N  =  1  KW  N  =0.1  Watt. 

X  3 


Na+Ni 


=  30  db.  The  interferer  carrier 


is  assumed  to  be  at  the  desired  signal  frequency.  If  the 
interferer  shifts  frequency,  his  angular  position  will  appear 
to  shift.  If  he  shifts  very  much,  part  of  his  signal  may  fall 
outside  of  the  signal  bandwidth. 

In  Fig.  9  a  second  30  db  Nj  interferer  is  added  at 
+120°  to  the  30  db  Nj  jammer  shown  in  Fig.  8.  a  second  minimum 
appears  as  expected.  These  interferers  are  far  enough  apart 
so  that  for  most  angles  of  signal  the  received  noise  level  is 
essentially  the  ambient  level,  and  (S/N)q  =  +16  db.  A  third 
interferer  could  be  added,  say  at  0°  without  changing  the  S/N)o 
for  the  existing  lobes  very  much.  However,  if  2  or  more  jammers 
are  close  together,  a  wider  combined  minimum  is  formed,  as  in 
Fig.  10. 

In  Fig.  10  three  Njyg  =  30  db  interferers  are  at  72°, 

90°,  and  108°.  The  ambient  ^  is  still  +10  db  at  each  filter  or 

wa 

+16  db  at  the  array  output.  This  could  only  represent  3  separate 
interferers,  each  at  the  given  angle  with  each  interfering  trans 
mitter  using  a  different  random  PN  sequence  with  the  carriers  all 
at  the  signal  carrier  frequency.  It  could  also  represent  3 
Gaussian  spectrum-matched  interferers. 


SNR  versus  angle  when  steering  vector  equals  signal  angle,  with  2  Inter ferers 
each  with  NT  -  103W,  one  at  20°,  one  at  45°;  N  »#1W,  S  «  1W. 


with  N_  *  10  W,  at  angles  72 
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In  Fig.  11  we  again  have  3  interferers  close  to 
each  other,  but  they  are  not  symmetric  about  90°.  Note  that 
the  minimum  SNR  values  are  not  affected  even  in  this  case, 

(all  -30  db) ,  but  the  range  of  angles  where  the  received  SNR  is 
within  5  db  of  the  ambient,  (+16  db) ,  is  much  reduced. 

Also  in  this  case  it  is  theoretically  possible  for 
a  single  interferer  at  one  of  the  3  angles  to  transmit  3 
frequencies  simultaneously,  with  the  frequencies  chosen  so  that 
the  phase  shifts  between  adjacent  elements  match  those  for  the 
assumed  geometry  shown  for  3  interferers.  The  adaptive  algorithm 
described  cannot  distinguish  between  these  two  cases.  The 
amount  of  shift  required  in  Fig.  11  would  probably  move  the 
shifted  frequencies  out  of  practical  receiver  bandwidths. 

But  smaller  shifts  must  be  considered  as  a  possible  jamming 
threat. 

Note  that  if  the  interferer  were  frequency  hopping 
between  the  3  frequencies  a  covariance  matrix  based  on  obser¬ 
vations  over  an  extended  time  interval,  as  required  by  some 
adaptive  narrow  band  algorithms,  would  contain  all  three 
frequency  contributions.  However,  our  measurements,  where  the 
covariance  matrix  is  based  on  wideband  observations  taken  within 
a  short-time  interval,  might  see  one  frequency  at  a  time  and 
shift  the  array  weights  to  track  the  apparent  change  in  inter¬ 
ferer  position. 

In  Fig.  12  the  performance  of  the  4  element  array 
changes  completely  as  one  increases  the  number  of  jammers  from 
3  interferers  to  4  interferers. 


Ttq,  11.  SNR  versus  angle  when  steering  vector  equals  signal  angle,  with  3  Interferers,  each 
with  N_  -  10^W,  at  angles  27°,  45°,  63°;  N  »  .1W;  S  «  1W. 
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Fig.  12.  SNBo  versus  angle  for  4-element  array.  Steering  vec  or 
9  equals  signal  angle  (top)  3  interferers,  with  the  N.  varied 

from  1,  10,  100,  to  1000W;  Na  =  1W;  S  =  1W?  (bottom!  4  Inter 
ferers,  with  other  conditions  same  as  for  top. 
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In  the  top  of  the  figure  we  have  3  interferers  with 
equal  powers,  one  at  30°,  one  at  90°  and  one  at  150°.  As  the 
interferers  powers  are  all  increased  from  20  db  Njys  to  50  do 
NI/S'  *n  the  received  SNR  form  in  the  jammer  directions 

and  increase  in  depth  with  increased  interferer  power.  How¬ 
ever,  two  lobes  in  the  direction  of  broadside  +  35°  remain. 


where  received  SNR  is  close  to  ambient,  i.e.,  regions  of  good 
interferer  suppression. 

In  the  bottom  half  of  the  figure  we  have  4  inter¬ 
ferers  roughly  equally  spread  in  angle  over  180°.  Now  as  the 
power  of  all  interferers  is  increased,  the  available  SNR  for  all 
received  directions  is  decreased.  For  the  geometry  shown  the 
decrease  in  available  SNR  is  approximately  equal  to  the  inter¬ 
ferer  power. 

Such  a  decrease  in  performance  is  of  course  expected 
when  the  number  of  interferers  exceeds  the  "number  of  nulls 


which  can  be  formed."  The  resulting  performance  often  is  not 
reported.  A  real  system  will,  however,  encounter  many  sources. 
In  particular,  an  interferer  may  try  to  generate  a  wide  band 
signal  such  that  more  than  one  "available  null"  is  used  up 
in  suppressing  a  single  jammer  location.  The  performance  under 
these  conditions  therefore  is  important  system  design  informa¬ 
tion. 


Fig.  13  to  18  explore  the  effect  of  a  4th  interferer 
in  more  detail.  In  all  these  cases,  3  of  the  interferers  are 
1  KW,  30  db  Nt/c  jammers.  For  Fig.  13  to  Fig.  17,  these  are 
located  at  45°,  90°  and  135°.  A  fourth  interferer  initially 
at  0°  is  increased  from  -10  db  Nj^s  to  30  db  Nj^s  as  we  go 
from  Fig.  13  to  Fig.  17.  The  shape  of  the  pattern  does  not 


Fig.  15.  SNRq  versus  angle  when  steering  vector  eguals  signal  angle,,  and  4  Interferers, 

3  of  which  have  NT  -  10^W  at  angles  45°,  90°,  and  135°;  4^  Interferer  at  0°  has  10W 


Fig.  16.  SNR  versus  angle  when  steering  vector  equals  signal  angle/  and  4  Interferers, 

3  of  which  have  N-  -  103W  at  angles  45°,  90°,  and  135°?  4th  Interferer  at  0°  has  100W 
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change  until  the  last  step,  but  it  is  clear  that  even  a  rather 
weak  4th  interferer  can  effectively  reduce  spatial  processing 
gain  for  all  angles. 

In  Fig.  18  a  weaker (lOO  T7 )  jammer  is  placed  at 
135°,  instead  of  at  0°  as  in  Fig.  17.  The  available  remaining 
pattern  is  quite  different  for  the  two  cases. 

These  figures  illustrate  some  of  the  factors  which 
must  be  considered  in  any  adaptive  array  design.  In  particular, 
it  is  important  to  consider  that  theoretically  there  are  wide 
range  of  angles  where  3  or  less  interferers  can  be  effectively 
suppressed  so  as  to  have  little  effect,  independent  of  inter¬ 
ferer  strength.  More  interferer  power  just  results  in  more 
effective  suppression.  This  communication  advantage  is  lost 
completely  if  one  or  more  additional  real  or  simulated  inter¬ 
ferers  are  introduced  into  the  covariance  matrix  noise  samples. 

The  system  must  be  designed  so  that  it  is  difficult  for  the  jammer 
to  produce  or  simulate  these  extra  interferers.  The  calculation 
in  these  sections  are  again  optimal  array  performance,  with 
assumed  complete  coherence  between  the  4  antenna  elements  for 
all  signals. 

We  must  emphasize  that  performance  associated  with  the 
above  curves  implied  a  large  number  of  apriori  steering  vectors 
or  beams  (one  for  each  angle  evaluated) .  While  in  concept  one 
could  envisage  a  large  number  of  apriori  beau*.  (or  steering 
vectors) ,  in  practice  we  are  as  usual  only  interested  in  the 
number  sufficient  to  give  the  bulk  of  the  achievable  gain. 

We  recall  that,  in  the  acquisition  mode,  both  time  of  signal 
arrival  and  direction  of  arrival  are  unknown.  We  also  recall 
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from  Sec.  3  that  in  the  case  of  the  time  function  we  used  a 
sliding  window  which  moved  in  increments  of  h  chip  duration 
and  detected  the  signal  if  it  exceeded  a  threshold,  at  one  or 
more  sliding  window  positions.  Similarly  we  can  choose  a  finite 
number  of  assumed  steering  vectors  for  each  time  window  and  de¬ 
tect  the  signal  when  the  combination  of  a  time  window  and  one  or 
more  steering  vectors  result  in  a  signal  output  above  a  threshold 
level.  To  determine  the  number  of  different  steering  vectors 
which  might  be  required  to  reliably  detect  the  signal  even 
when  it  does  not  match  the  assumed  steering  vectors  exactly,  - 
we  must  know  the  beam  width  for  different  steering  vectors 
and  different  numbers  of  interferers. 

Again,  we  do  not  have  an  equation  for  the  antenna 
pattern  when  an  arbitrary  number  of  interferers  are  present, 
but  we  do  have  computer  programs  for  calculating  the  antenna 
pattern  for  the  4  element  linear  array  and  any  assumed  inter¬ 
ferers  pattern.  The  examples  given  below  were  obtained  by 
the  CMAT  5  program  listed  in  Appendix  B.  The  antenna  array 
weights  are  normalized  so  that  the  noise  output  is  unity.  The 
decibel  value  plotted  therefore  match  the  performance  SNR  when 
the  signal  coincides  with  the  steering  direction. 

Starting  in  Fig.  19  we  first  show  a  number  of  ambient 
noise  antenna  patterns  which  are  generated  by  the  array  when  no 
directional  interferers  are  present.  In  Fig.  19  the  steering 
vector  is  broadside  at  +90°.  The  4  element  pattern  has  3  lobes 
with  nulls  as  shown.  We  consider  the  steering  vector  nulls  to 
be  perfect  nulls  (zero  gain),  where  the  phase  of  the  received 
signal  changes  phase  by  180°.  This  is  in  contrast  to  the  minima 
of  finite  depth  in  the  direction  of  a  interferer  which  may  be 
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produced  by  the  adaptive  algorithm.  For  our  antenna  array  model 
such  perfect  steering  vector  nulls  exist.  A  real  antenna  array 
may  have  defects  such  that  the  in-phase  and  quadrature  signals 
do  not  go  through  zero  at  the  same  angle.  It  is  important  to 
distinguish  between  these  nulls  of  the  steering  vector  pattern, 
and  "nulls"  formed  by  the  adaptive  algorithm.  The  patterns  in 
the  lower  half  plane  are  mirror  images  of  the  patterns  shown. 

Figs.  20  and  21  show  similar  ambient  beam  patterns  for 
steering  angles  of  45°  and  0°,  respectively.  The  135°  beam  pat¬ 
tern  would  be  a  mirror  image  of  the  45°  pattern.  Note  that  by 
using  only  2  steering  vectors,  45°  and  135°,  provides  good 
coverage  for  signal  reception  is  given  for  all  directions 
except  the  nulls  near  broadside.  Adding  a  third  steering 
vector  of  90°  covers  these  nulls.  We  could  also  add  a  4th 
steering  vector  at  0P  for  some  redundant  coverage.  The 
conclusions  are  that  about  4  steering  vectors  is  enough  to 
provide  360°  reception  in  ambient  noise. 

2 

Fig.  22  shows  the  antenna  power  pattern  *  when  a 
interferer  at  90°  is  added  to  the  45°  steering  vector  ambient 
pattern.  A  minimum  or  "null"  more  than  50  db  deep  forms  in 
the  interferer  direction.  The  rest  of  the  lobe  structure 
also  changes.  However,  the  45°  pattern  plus  the  135°  pattern 
still  provide  essentially  complete  360°  coverage  for  signal 
array  from  the  interferer  at  90°  (and  at  -90°) . 

Fig.  23  adds  another  30  db  interferer  at  60°,  while 

1.  in  most  practical  applications  a  linear  array  would  be  mounted 
on  a  ground  plane  which  suppresses  the  lower  half  space. 

The  received  noise  is  always  normalized  to  one;  therefore,  the 
antenna  gain  pattern,  which  specifies  the  arriving  signal  power 
is  equal  to  the  output  signal-to-noise  ratio. 
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Fig.  23.  SNRq  versus  angle  for  a  steering  vector  at  45  ,  with  two  Interferers  of  Nj  ■  10  W 
each,  at  60°  and  90°,  and  N  *  S  *  1W. 


steering  vector  at  45°;  with  three  Interferers  of 
90°,  and  150°,  and  N  -  S  «  1W. 
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Fig.  24  adds  a  third  30  db  interferer  at  150°.  There  are 

still  large  lobes  which  would  give  good  signal  reception  for 
signals  not  near  these  interferers. 

In  Figs.  25  and  26  we  have  added  two  and  three 
interferers  to  the  135°  steering  vector  pattern.  The  com¬ 
bination  of  this  with  the  corresponding  45°  patterns  gives  at 
least  one  pattern  with  good  coverage  for  all  angles  away  from 
the  jammers. 

Figs.  27,  28  and  29  show  something  different.  The 
steering  vector  here  is  0°.  Adding  the  two  30  db  inter¬ 

ferers  at  90°  and  60°  in  this  case  appears  to  have  no  effect  at 
all  on  the  0°  steering  angle  ambient  pattern.  The  ambient 
pattern  for  this  steering  vector  already  had  infinite  nulls  at 
90°  and  60°. 

Adding  the  3rd  jammer  does  introduce  a  new  minimum. 

In  spite  of  the  fact  that  these  jammers  are  in  "nulls"  of  the 
beam  pattern,  the  covariance  matrix  does  include  the  terms  to 
form  adaptive  nulls  in  these  directions.  If  a  4th  interferer 
is  introduced^'  the  response  degrades  just  as  badly  as  if  no 
begun  nulls  had  occurred  in  any  of  the  jammer  directions.  The 
algorithm  does  not  appear  to  be  exploiting  the  steering  vector 
nulls  efficiently  in  this  particular  case. 

Figs.  30,  31  and  32  show  this,  plus  another  effect. 
Here  we  happen  to  have  chosen  a  steering  vector  of  90°  which 
coincides  with  the  direction  of  the  90°  interferer.  The  90° 
jammer  in  this  case  can  not  produce  a  null,  instead  it  reduces 
the  antenna  gain  in  all  directions.  If  patterns  for  other 
steering  vectors  are  available,  this  may  not  be  serious. 
Otherwise  communication  for  all  directions  is  denied  by  the 


1.  No  Figure  is  included  to  portray  this  4-interferer  case. 


Fig.  25.  SNR0  versus  angle  for  a  steering  vector  at  ] 
Nt  »  103W  each,  at  60°  and  90°,  and  N  -  S  = 


Pig.  26.  SNRq  versus  angle  for  a  steering  vector  at  135°,  with  three  Interferers  of 
N_  *  103W  each,  at  60°,  90°,  and  150°,  and  N  -  S  *  1W. 
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steering  vector  at  0  P  with  one  Interferer  of 


Pig.  29.  SNRq  versus  angle  for  a  steering  vector  at  0°,  with  three  Interferers  each  with 
N_  «  103W  at  60°,  90°,  and  150°,  and  N  -  S  =  1W. 
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Pig.  32.  SNRq  versus  angle  for  a  steering  vector  at  90  ,  with  three  Interferers  each  with 
N,  -  103W  at  60°,  90°,  and  150°,  and  N  -  S  *  1W. 
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interferer  which  happened  to  be  located  in  the  chosen  (assumed) 
steering  vector  direction,  even  though  the  signal  may  not  be  in 
the  interference  direction. 

In  Fig.  31  the  60°  jammer  is  added.  This  happens  to 
fall  in  a  steering  vector  null  and  does  not  further  degrade 
the  pattern.  The  150°  jammer,  when  added,  (Fig.  32)  does  change 
the  pattern. 

We  see  from  the  last  example,  that  in  choosing  a 
minimum  set  of  steering  vectors  one  should  have  some  redundancy. 
When  one  of  the  steering  vectors  is  essential  for  coverage  of 
certain  signal  arrival  directions,  an  interferer  which  becomes 
too  near  that  steering  vector  can  reduce  the  coverage  from  that 
steering  vector  and  leave  a  gap  where  no  signal  can  be  received. 
This  gap  is  in  addition  to  the  gap  which  always  forms  near  the 
jammer  direction. 

Figs.  33,  34,  35  and  36  show  that  the  algorithm,  in 
general,  does  not  reduce  the  antenna  power  gain  to  zero  in  the 
direction  of  the  interferer,  but  only  reduces  the  gain  by  some 
rather  large  amount  in  the  jammer  direction.  The  reduction  in 
gain  is  chosen  so  as  to  maximize  the  output  SNR.  When  the  inter¬ 
ferer  is  much  larger  than  either  ■►he  signal  or  the  ambient  noise, 
then  the  optimum  process  requires  a  deep  null  at  the  jammer.  If 
the  interferer  is  weak,  other  factors  become  more  important  for 
maximizing  output  SNR.  Note  again  in  Fig.  35  that  if  by  accident 
an  interferer  happens  to  line  up  with  an  assumed  steering  vector, 
the  gain  for  that  steering  vector  is  suppressed  in  all  directions 


Fig.  34.  SNRo  versus  angle  for  a  steering  vector  at  45  ,  with  three  Interferers  of  Nj  =  10W 
each  at  60°,  90°,  and  150°,  and  N  »  S  »  1W. 
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Figs.  37,  38,  39  and  40  were  generated  to  determine 

how  close  to  the  steering  angle  a  jammer  must  be  in  order  for 

the  interferer  to  completely  sup;-  =ss  the  antenna  pattern  in 
all  directions.  It  appears  that  for  suppression  in  all  directions 

the  interferer  must  be  within  a  degree  or  so  of  the  chosen  steer¬ 

ing  angle  for  our  4  element  array.  Otherwise  only  a  section  of 
the  ambient  pattern  is  suppressed.  Given  4  fixed  steering 
angles,  it  is  unlikely  that  more  than  one  of  3  jammers  would 
fall  on  any  of  the  steering  angles.  However,  there  often  are 
many  other  low  level  sources  in  the  antenna  field.  All  of  these 
potentially  can  reduce  the  adaptive  array  gain,  but  only  in 
proportion  to  their  strengths. 

Further  examples  demonstrate  that  5  steering  vectors, 
leads  to  acceptable  losses  at  least  during  acquisition  from  the 
ideal  performance  which  would  be  possible  if  the  steering  vector 
were  directed  at  the  actual  signal  source. 

Fig.  41  assumes  that  the  weights  have  been  calculated 
for  4  steering  directions,  0,  45,  90,  and  135°.  The  signal  has 
been  processed  for  each  of  the  weights.  If  the  signal  to  noise 
output  exceeds  a  preset  threshold  on  any  of  the  processed  outputs, 
a  signal  is  considered  to  be  detected.  The  signal  is  then  de¬ 
modulated  to  see  if  a  useful  message  is  present.  In  Fig.  41, 
the  solid  line  shows  the  improvement  in  SNR  due  to  the  array 
as  a  function  of  signal  angle  when  no  interferers  are  present. 

(The  larger  gains  of  the  4  steering  directions  (dotted  lines)  is 
always  chosen) .  Close  to  the  expected  6  db  gain  of  the  array 
is  available  over  most  angles.  There  is  loss  of  =  2db  over 
omnidirectional  coverage  near  20°  +  110°.  This  would  be  reduced 
by  changing  the  45°  +  135°  steering  angles  to  angles  closer  to 
broadside.  The  ideal  curve  for  steering  angles  always  matched 


Pig.  37.  SNRq  versus  angle  for  a  steering  vector  of  85°,  with  three  Interferers  of  Nj  *  10W 
each  at  60°,  90°,  and  150°,  and  N  *  S  »  1W. 


Fig.  39.  SNRq  versus  angle 
Nt  -  10W  each  at  6 
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to  signal  direction  would  be  a  6  db  circle. 

When  three  interferers  are  present,  the  SNR  available 
with  these  4  steering  vectors  is  shown  in  Fig.  42.  Note  that 
in  this  case  one  of  the  interferers  happened  to  coincide  in 
direction  with  the  90°  steering  vector.  This  resulted  in  a 
large  reduction  in  the  antenna  gain  available  for  this  steering 
vector  (shown  in  dashed  lines) .  If  the  90°  interferer  had  been 
as  much  as  2°  or  5°  away  from  the  90°  steering  vector,  the 
antenna  gain  for  the  90°  steering  vector  would  have  been  much 
higher  for  all  signals  except  those  coming  from  the  interferer 
direction.  Fig.  43  shows  the  performance,  for  the  same  inter¬ 
ferer  conditions  as  Fig.  42,  except  that  the  correct  steering 
vector  is  used  for  every  angle. 

We  see  that  the  worst  performance  with  4  preset 
steering  factors  occurs  when  an  interferer  happens  to  be  exactly 
aligned  with  one  of  the  preset  steering  angles.  If  the  signal 
also  happens  to  arrive  from  an  angle  near  the  interferer,  its 
SNR  is,  of  course,  reduced  both  in  the  ideal  case  and  when 
using  only  4  steering  vectors.  However,  with  only  4  steering 
vectors,  the  SNR  for  some  signals  not  near  the  interferer  are 
also  reduced  somewhat  from  the  ideal  curve.  This  occurs  only 
if  the  interferer  is  within  +  2°  to  +  5°  of  one  of  the  steering 
angles,  which  occurs  with  low  probability.  To  avoid  this 
potential  (minor)  problem,  we  suggest  using  5  preset  steering 
angles.  If  an  interferer  falls  close  to  one  of  these,  there 
will  still  be  4  other  patterns  which,  combined,  give  '-lose  to 
theoretical  optimum  performance. 


versus  angle  when  the  best  of  steering  vectors,  at  0  ,  45  ,  90  ,  and  135 
ised,  with  three  10^W  Interferers  at  60°,  90°,  and  150°,  and  N  =  S  =  1W. 


SNRq  versus  angle  for  the  same  Interferer  conditions  as  in  Pig.  42,  where  the 
correct  steering  vector  is  now  used  for  each  angle. 
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The  conclusions  are: 

(1)  The  algorithm  produces  antenna  patterns  (equivalent  to 
SNR  patterns)  which  in  the  absence  of  interferers  give 
quite  broad  receiving  sectors.  Three  or  four  steering 
vectors  cover  360°  adequately. 

(2)  The  presence  of  interferers  introduces  narrow  minima/ 
not  nulls,  in  the  ambient  patterns  in  the  direction  of 
the  interferers.  The  depth  of  the  minimum  Is  automati¬ 
cally  adjusted  to  maximize  the  received  SNR  when  the  signal 
comes  from  the  direction  of  the  steering  vector. 

(3)  The  received  SNR  is  good  over  wide  angles,  both  near  the 
steering  angle  and  in  a"sidelobe  area." 

(4)  Predicting  what  will  happen  in  the  real  world  is  compli¬ 
cated  because  (a)  interferers  can  appear  in  nulls  of 
patterns  formed  by  other  interferers;  (b)  interferers  can 
appear  at  the  prechosen  steering  angles.  In  the  first 
case  the  new  interferer  has  no  effect  in  pattern  on  SNR, 
but  does  use  up  a  potential  null  produced  by  the  co- 
variance  matrix.  In  the  second  case,  the  pattern  in  all 
directions  is  suppressed  for  that  steering  vector. 

4.3  Effect  of  Number  of  Samples  Used  In  Estimating  Matrix 
The  results  in  Section  4.1  and  Section  4.2  are 
theoretical  and  limiting  in  the  sense  that  it  is  assumed  the 
correct  covariance  matrix  is  known;  i.e.,  the  spatial  angles 
and  interferer  and  signal  powers  were  inserted  as  analysis 
input.  In  the  operational  case,  the  covariance  matrix  must 
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be  estimated  from  EM  signals  arriving  at  the  elements.  A 
basic  item  of  interest  is  how  many  independent  sets  of  samples 
of  array  output  are  needed  to  adequately  estimate  the  covari¬ 
ance  matrix.  We  first  consider  the  ideal  case  where  the  steer¬ 
ing  vector  is  aligned  with  the  desired  signal  direction. 

The  analysis  and  conclusions  in  Ref.  9  directly  address 
this  issue.  The  array  processing  envisaged  in  Ref.  9  consisted 
of  both  time  and  space  sampling  for  SMI  directly  at  the  antenna 
elements.  However/  with  intervening  matched  filters,  as  in  our 
case,  the  filter-output  samples  retain  their  spatial  integrity, 
even  though  the  individual  samples  themselves  are  transferred 
by  the  matched-filter  transfer-function.  Therefore  one  should 
expect  that  the  results  concerning  number  of  samples  needed  to 
estimate  the  covariance  matrix  reached  in  Ref.  9,  should  apply 
to  our  case. 

We  will  briefly  restate  the  results  here.  Ref.  9 
uses  a  highly  theoretical  development  to  show  that  the 
normalized  signal-to-noise  ratio  (rates  of  output  SNR  with 
ideal  covariance  matrix  to  output  SNR  with  estimated  matrix 
behaves  as: 

Normalized  SNR  (db)  =  10  log^g  '^^+1'^"'  (18) 

where:  n  -  number  of  independent  samples  used  in  esti¬ 
mating  covariance  matrix. 

K  =  number  of  elements  in  the  array 
If  one  examines  the  3  db  loss  point,  one  finds  n  =2K-3.  Hence, 
n ~  2K  if  K  >  10,  for  3  db  loss. 

We  should  note  that  Equation  (18)  implies  that  only 


one  sample  is  taken  per  element  per  sampling  epoch;  in  the 


general  case,  without  matched  filters  ahead  of  the  SMI  process¬ 
ing,  one  will  wish  to  take  (say)  Z  samples  per  element,  each 
separated  by  a  small  delay  compared  to  the  independent  sample 
time,  in  order  to  prevent  decorrelation  of  samples  from  wide¬ 
band  interferers.  In  that  case,  the  n  would  be  n  =  2k£.-3  £2  2kJl. 

In  our  case,  with  matched  filters  ahead  of  the  SMI  processing, 
there  is  much  less  advantage  in  taking  more  than  one  sample  per 
element  because,  at  current- (matched)  sample  time,  the  band¬ 
width  is  collapsed  to  that  of  the  information  bandwidth  (see 
Section  5.2). 

The  result  in  Equation  (18)  was  derived  by  finding  the 
probability  distribution  of  the  normalized  SNR  (using  a  complex 
Wishart  distribution) ,  and  verifying  the  result  using  simula¬ 
tion  where  30  scatterers  (as  in  radius  clutter)  provide  the 
interference.  In  reviewing  these  results,  our  interest  centered 
in  the  small  element  (K  =  4)  and  small  number  of  interferer 
(nj  5  3)  case. 

First  we  should  note  that  the  relatively  small  number 
of  required  samples  from  Equation  (18)  is  surprising  in  the 

A 

sense  that,  with  these  sample-sizes,  the  actual  M  is  not  close 
to  M  because  of  the  cross-products  resulting  from  the  ambient 
and  the  large  interferer  terms.  Upon  analyzing  this  situation, 
we  could  show  (see  Ref.  13)  that  the  important  requirement  is 

A 

that  the  M  eigenvalues  line  up  with  M  eigenvalues.  Therefore, 
only  a  relatively  small  number  of  samples  may  be  required. 

Fig.  44  shows  degradation  in  normalized  SNR  as  a 
function  of  number  of  samples  n,  and  number  of  antenna  elements  K. 
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It  is  based  on  Equation  (18)  .  For  4  element  and  8  samples, 
the  expected  loss  is  1.6  db.  The  number  of  samples  required 
increases  rapidly  if  better  performances  are  required,  or 
larger  numbers  of  elements  are  used..  In  Section  6,  we  report 
the  results  of  a  computer  program  which  does  both  the 
matched  filter  time  processing  and  the  SMI  array  processing. 

The  examples  which  have  been  computed  appear  consistent  with 
the  above  theory.  More  examples  would  be  needed  for  a  statisti¬ 
cal  proof. 

When  the  steering  vector  is  not  aligned  with  the 
received  signal,  the  analysis  is  more  complicated  (Ref.  11) . 

In  general,  more  samples  are  needed  for  a  given  level  of  per¬ 
formance.  The  number  of  samples  required  depends  on  the  de¬ 
gradation  introduced  by  using  the  approximate  steering  vector. 
Since  the  degradation  is  small  in  our  examples,  we  do  not  ex¬ 
pect  a  significant  increase  in  the  number  of  samples  required. 
However,  we  have  no  data  to  confirm  this. 
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5.  Comparative  Advantage  of  PN  Matched  Filter  Plus  SMI 

We  will  compare  the  performance  here  (using  Fig.  6) 
with  the  general  performance  of  the  LMS  closed-loop  method 
depicted  in  Fig.  3  and  the  straight-forward  method  of  using 
a  power  inverting  adaptive  array  followed  by  a  matched  filter. 
These  three  appear  to  be  the  clear  logical  choices  for 
combining  PN  with  adaptive  array. 

In  a  communications  application  involving  spread 
spectrum  and  adaptive  arrays,  the  system  must  first  acquire 
synchronization  in  both  time  (for  the  temporal  processing) 
and  direction  (for  the  spatial  processing)  before  the  expected 
resistance  to  interference  occurs.  Therefore,  it  is  the 
behavior  in  the  acquisition  mode  that  is  the  driving  limitation 
in  any  system  that  combines  spread  spectrum  with  adaptive 
processing.  First  we  recall  that  the  asymptotic  steady-state 
response  of  an  LMS  algorithm  with  known  reference  uses  weights 
which  are  proportional  to  the  same  inverted  noise  covariance 
matrix  that  is  estimated  and  used  directly  in  any  SMI  processing. 
This  leads  us  to  conclude  that  the  LMS  loop  of  Fig.  3  and  our 
procedure  of  Fig.  6  should  give  the  same  nominal  steady-state 
results  if:  1)  the  PN  modem  in  the  reference  loop  is  syn¬ 
chronized;  and  2)  only  narrowband  or  CW  interference  is  en¬ 
countered.  We  will  rgue  that  the  PN  matched  filter  arrange- 
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ment  using  multiple  steering  vectors  is  superior  in  the 
acquisition  mode  (matched  time  and  direction  unknown)  and 
when  matched- spectrum  interference  is  contemplated. 

First  consider  the  case  where  an  adaptive  array 
(as  in  Fig.  2)  is  followed  by  a  PN  matched  filter.  In  this 
case  the  adaptive  array  cannot  separate  signal  samples  from 
noise  samples  (for  a  continuous  signal) ,  and  the  array  will 
treat  the  intended  signal  as  another  interferer:  a  signal 
will  be  nulled  in  accordance  with  its  strength.  Such  "power 
inversion"  algorithms  (Ref.  10)  can  succeed  if  the  ensueing 
temporal  processing,  via  the  matched  filter,  can  restore  a 
nulled  signal  to  an  adequate  minimum  output  SNR.  Even  if 
this  occurs,  nulling  the  signal  has  two  significant  dis¬ 
advantages  over  either  the  Fig.  3  closed- loop  method  or  our 
method  (Fig.  6) :  1)  degrees  of  freedom  are  consumed  by  the 
nearby  (intended)  signal  transmitters,  anu  2)  the  field-of- 
view  of  the  receiver  has  a  segment  or  slice  removed,  that  may 
deny  access  to  another  signal  that  is  in  or  near  the  null. 

Therefore,  we  believe  it  is  basic  to  attempt  to 
prevent  nulling  the  signal  when  combining  PH  spread  spectrum 
with  adaptive  arrays.  Since  the  methods  of  Fig.  3  and  Fig.  6 
both  have  this  property,  we  will  emphasize  their  comparative 
properties. 

It  is  necessary  to  emphasize  that  the  significant 
difference  between  the  two  systems  does  not  focus  during 
steady- state,  but  occurs  in  the  manner  and  success  with  which 
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they  acquire  in  the  acquisition  mode,  under  worst  case  con¬ 
ditions.  We  will  compare  the  two  approaches  in  two  aspects: 
acquisition  and  bandwidth  effects. 

5.1  Acquisition  Comparison 

Consider  first  the  spatial  uncertainty  in  the 
acquisition  mode.  Using  a  matched  filter  followed  by  SMI 
affords  a  natural  way  to  handle  spatial  uncertainty:  use 
parallel  processing  consisting  of  multiple  apriori  steering 
vectors  corresponding  to  hypothesized  signal  directions, 
utilizing  digital  processing.  In  the  above  section  we  demon¬ 
strated  by  a  series  of  experiments  that  one  needs  only  4  or  5 
such  apriori  beams  if  the  number  of  elements  is  4. 

In  effect,  closed-loop  algorithms,  such  as  in  Fig.  3, 
tend  to  naturally  be  "serial"  and  analog,  while  open-loop  ones 
extend  naturally  to  parallelism  and  open-ended  digital  processing. 

The  next  pertinent  issue  is  transient  response  of 
the  adaptive  processing.  As  is  well-known,  the  time-constant 
of  closed-loop  schemes  depends  on  both  the  loop  gain  and  the 
level  of  the  received  interfering  power.  T1 '  response  time  is 
limited  by  filter  bandwidths  and  loop  stability  requirements. 

As  a  result,  the  performance  of  a  scheme  as  in  Fig.  3  is  dependent 
on  the  actual  signal,  ambient,  and  interference  levels  encountered. 
With  our  use  of  the  SMI  method,  convergence  is  independent  of  the 
interference  environment. 
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Whereas  conventional  SMI  will  require  use  of  2K-3 
consecutive  independent  (waveform)  samples,  note  that  with 
matched  filters  one  can  compute  the  necessary  samples  used 
for  estimating  the  covariance  matrix  for  the  weights  at  each 
current  sample,  using  the  current  contents  of  the  transversal 
filter,  by  using  a  rapid  digital  processor.  This  is  only 
slightly  different  from  using  the  (say)  eight  previous  filter- 
output  samples  to  evaluate  the  weights  for  the  current  sample. 

In  any  event,  the  transient  response  for  either 
approach  is  essentially  instantaneous,  and  independent  of  the 
received  levels.  This  should  be  especially  important  in 
practice,  especially  for  interferers  that  would  attempt  special 
methods  such  as  blinking. 

We  now  consider  temporal  synchronizing.  To  be  truly 
interference-resistant,  the  PN  modem  in  Fig.  3  must  utilize 
a  "long-code",  so  that  time-sync  must  be  searched  over  some 
apriori-established  time-uncertainty.  As  noted  in  Ref.  2, 
the  traditional  way  to  accomplish  time-sync  is  to  use  code 
slewing,  where  the  local  PN-code  generator  is  run  faster 
(or  slower)  than  the  incoming  signal  code.  Hence,  relative 
to  the  onset  of  intended  signal  transmission,  there  is  added 
delay  incurred  between  signal-onset  and  code-alignment.  This 
causes  what  is  sometimes  called  the  "great  race"  problem:  if 
the  array  time  constant  is  faster  than  the  modem  acquire  time- 
constant,  the  array  will  null  the  signal  before  it  can  be 
acquired  by  the  modem. 
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This  delay  does  not  occur  with  our  arrangements,  be¬ 
cause  the  matched-filter  approach  implies  that  a  certain  PN 
segment  is  anticipated,  and  is  set-up  in  the  filter.  The 
matched-filter  must  anticipate  maximum  range-uncertainty  plus 
clock  drift.  Currently  emerging  CCD  and  SAW  programmable 
matched  filters  permit  examing  different  range  bins  at  an 
internal  processing  speed  that  is  higher  than  the  spread- 
spectrum  bandwidth.  In  any  event,  so  long  as  some  method  is 
found  to  provide  two^*  full- segments  of  a  PN  sequence,  there 
will  be  no  additional  delay  due  to  time  sync,  after  onset  of 
signal.  One  can  avoid  interference-exploitation  of  the  two- 
period  transmission  by  making  it  impractical  to  predict  the  PN 
segment  that  will  be  used. 

5.2  Performance  With  Wideband  Interference 

A  key  issue  that  must  be  emphasized,  when  combining 
adaptive  arrays  with  spread  spectrum,  concerns  the  degraded 
array  performance  with  the  conventional  single-tap  approach 
(as  in  Fig.  3)  when  the  interferer  uses  a  spectrum-matched 
interference  signal.  As  one  increases  the  intended  signal 
bandwidth  via  spread-spectrum  to  increase  temporal  processing 
gain,  one  of  course  must  anticipate  possible  use  of  an  inter¬ 
ference  signal  of  equal  bandwidth. 

When  array  processing  precedes  temporal  processing, 
as  in  Fig.  2,  a  de-correlation  of  interference  signals  occurs 
due  to  the  finite-length  of  the  array,  at  all  incidence  angles 
except  broadside,  if  one  uses  only  one  tap  per  antenna  element. 

1.  It  is  possible  to  avoid  requiring  two  identical  segments, 
but  this  is  beyond  the  scope  of  the  effort  here. 


This  will  cause  the  array  processing  gain  to  decrease  linearly 
with  increase  in  (PN)  signal  bandwidth.  Adding  multiple  taps 
to  each  element  as  described  in  Ref.  3  and  Ref.  5  can  essen¬ 
tially  remove  the  de-correlation  effect. 

One  can  demonstrate  the  result  of  de-correlation, 
in  the  steady-state,  by  assuming  an  interference  spectral 
shape,  finding  its  autocorrelation,  and  then  computing  the 
expected  value  of  array  cancellation  as  a  function  of  inter¬ 
ference  bandwidth  and  maximum  array  time-delay  (for  end-fire, 
if  a  linear  array)  ,  as  was  done  in  Ref.  5.  In  Appendix  E 
we  calculated  the  de-correlation  effect  if  only  a  single  tap 
is  used  by  first  assuming  an  autocorrelation  function  for 
the  interference  and  a  single  interferer.  Then,  using  the 
(presumed  known)  corresponding  covariance  matrix  of  the  inter¬ 
ferer,  the  smallest  eigenvalue  (which  is  the  steady-state 
solution)  is  found,  as  a  function  of  bandwidth.  In  general, 
it  was  found  that,  with  the  single  tap",  the  array  processing  gain 
decreases  linearly  with  increase  in  interference  (and  PN  signal) 
bandwidth  (see  Appendix  E) .  As  noted,  the  chief  consequence  of 
this  issue  is  that  the  system  of  Fig.  3  must  be  altered  to  in¬ 
clude  multiple  taps  per  antenna  element. 

To  see  how  our  approach  (see  Fig.  6)  or  a  "code¬ 
stripping"  approach,  as  in  Fig.  4,  would  fare  in  this  de-correlation 
aspect,  we  assumed  one  interferer  (see  Appendix  E)  and  a  broadside 
signal  with  unity  power.  The  process  output,  3(t)  was  then 
formed,  in  terms  of  the  output  interference  power.  Similar  to 
the  procedure  used  in  Ref.  5,  we  assumed  an  autocorrelation 
formation  for  the  filter  output  samples  due  to  interference. 
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We  then  assumed  the  interference  bandwidth  matches 
that  of  the  signal,  and  evaluated  the  minimum  eigenvalues, 
as  before.  We  found  that,  using  the  assumed  exponential 
autocorrelation,  that  the  interference  output  power  is 
significantly  lower  than  was  the  case  for  the  previous 
(single  tap)  conventional  arrangement  (Fig.  3.).  Al¬ 
though  this  was  done  for  a  single  interferer,  we  believe  similar 
results  will  hold  for  multiple  interferers,  up  to  the  permissible 
degrees  of  freedom.  Furthermore,  the  advantage  for  temporal 
processinc  ahead  of  spatial  processing  shown  in  Ref.  14  is 
dependent  on  the  nature  of  the  RF  filter  that  could  precede  the 
temporal  processing. 

We  conclude,  then,  that  if  the  closed-loop  arrange¬ 
ment  used  only  a  single  tap  per  element,  as  in  Fig.  3,  that 
our  arrangement  (as  in  Fig.  6)  would  outperform  it  by  an  amount 

I 

dependent  on  filter  assumptions.  However,  if  the  closed-loop 
arrangement  included  multiple  taps,  it  is  likely  that  the 
steady-state  performance  of  the  two  arrangements  would  be  essen¬ 
tially  equal. 
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6 .  Simulation  Demonstration 

We  constructed  a  computer-based  demonstration  of  the 
arrangement  shown  in  Fig.  6  both  to  provide  evidence  that  the 
process  works,  and  to  serve  as  a  learning  aid.  Now  we  esti¬ 
mate  the  covariance  matrix  from  simulated  noise  samples, 
rather  than  assume  it  is  known  (as  in  Sec.  4.2).  A  PN  matched- 
filter  of  length  255  was  used  for  the  temporal  processing,  and 
a  4-element  linear  array  spaced  at  j  was  used.  Six  apriori 
steering  vectors  were  used,  placed  at  angles  of  +  15°  +  45°  and 
+75°  from  broadside.  The  adaptive  array  output  for  each  of 
these  was  computed  on  a  sliding  window  basis.  The  outputs  are 
compared,  with  a  threshold,  looking  for  signal  peaks.  An  es¬ 
timated  output  SNR  (based  on  a  limited  number  of  noise  samples) 
was  also  calculated  for  each  steering  vector. 

The  desired  signal  consisted  of  2  or  more  periods 
of  a  255  bit  maximal  sequence.  This  signal  was  assumed  to 
phase  modulate  a  UHF  carrier.  The  source  was  assumed  to  be 
located  at  +15°,  which  corresponds  to  one  of  the  apriori  steering 
vectors.  The  phase  of  the  received  signal  at  each  antenna  element 
was  calculated.  The  computer  was  used  to  generate  the  complex 
samples  of  the  demodulated  signal  for  each  antenna  at  the  rate 
of  4  complex  samples  per  sequence  bit.  Similar  complex  samples 
for  the  interference  and  ambient  noise  were  added  to  the  signal 
samples. 

The  matched  filter  outputs  for  filters  consisting  of 
one  period  of  the  sequence  were  calculated  for  the  in-phase  and 
quadrature  component  at  each  antenna  element  on  a  sliding  window 
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basis.1'  The  window  moves  in  increments  of  a  sequence  bit  to 
avoid  missing  a  signal  peak.  The  resulting  set  of  complex  numbers 
are  the  signal  inputs  to  the  adaptive  array  processing. 

One  or  two  interfering  sources  were  simulated  by 
random  number  generators.  They  were  assumed  to  transmit  car¬ 
riers  phase  and  amplitude  modulated  by  Gaussian  noise.  For 
assumed  interference,  source  powers  and  source  positions,  the 
received  amplitudes  and  phases  were  calculated.  Corresponding 
random  in-phase  and  quadrature  demodulated  interference  samples 
were  generated  in  the  computer  and  added  to  the  corresponding 
received  signal  samples.  For  each  antenna  element,  independent 
in-phase  and  quadra tive  samples  of  Gaussian  noise,  representing 
ambient  noise,  were  also  generated  and  added  to  the  signal  samples. 

The  signal-free  complex  noise  samples  for  calculating 
each  sample  covariance  matrix  were  obtained  from  the  outputs  of 
the  matched  filter  with  the  window  displaced  from  the  correct 
matched  signal  output.  From  1  to  200  samples  of  the  covariance 
matrix  were  calculated  and  averaged  to  investigate  the  effect  of 
the  number  of  samples  used  in  the  estimate  of  the  covariance  matrix. 
The  nearest  noise  sample  was  displaced  2  bits  from  the  expected 
signal  output  and  the  additional  samples  as  required  were  dis¬ 
placed  further  in  increments  of  h  bit;  i.e.  8  signal-free  noise 
samples  were  obtained  by  sampling  between  2  bits  and  6  bits  from 
the  expected  signal  position  at  2  complex  samples  per  bit.  For 

1.  The  autocorrelation  function  of  a  maximal  sequence  of  length  L 
has  a  peak  of  amplitude  L  but  is  not  exactly  zero  between 
peaks.  It  can  be  made  exactly  zero  by  using  a  phase  modulation 
slightly  less  than  180°.  It  can  also  be  made  exactly  zero  by 
using  a  modified  "matched  filter":  a  "replica  reference"  con¬ 
sisting  of  only  the  L-l  "ones"  in  the  sequence.  This  was  done 

2 

in  the  program  described  here. 
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each  noise  sample  position,  a  sample  covariance  matrix  was 
calculated  between  the  samples  from  the  different  antenna 
elements . 

The  estimated  covariance  matrix  was  obtained  by 
averaging  the  sample  matrices.  It  was  then  inverted  and 
multiplied  by  one  of  the  apriari  steering  vectors  to  find 
one  set  of  combining  weights.  The  matched  filter  signal 
output  was  combined  using  these  weights  as  one  possible 
signal  output.  This  process  was  repeated  for  each  of  the  6 
assumed  steering  vectors.  The  complete  process  from  finding 
a  new  covariance  matrix  to  calculating  6  more  possible  signal 
outputs  was  repeated  continuously  on  a  sliding  window  basis. 

While  this  general  program  is  flexible  and  can  be 
used  for  a  number  of  exploratory  purposes,  we  will  here  empha¬ 
size  three  basic  aspects:  1)  the  behavior  of  the  array-pattern 
as  a  function  of  the  number  of  samples  used  in  estimating  the  co- 
variance  matrix,  for  various  assumed  conditions  and  2)  the  be¬ 
havior  of  the  output  time-axis  in  the  vicinity  of  the  peak  (for 
various  conditions)  to  demonstrate  typical  behavior. 


6. 1  Behavior  of  Array-Patterns  Versus  N 

Here  we  report  how  the  simulation  array  patterns 
behave,  as  a  function  of  the  number  of  samples  used  in  esti¬ 
mating  the  covariance  matrix.  The  series  of  results  portrayed 
in  Fig.  45  through  Fig.  49  describe  the  array  pattern,  which 
is  again  equivalent  to  an  output  SNR,  for  one  interferer 
through  5  interferers,  respectively.  For  each  Figure,  the 
results  are  shown  for  estimating  the  covariance  matrix  based 
on  either  n  =  2,4,8,16,32,200.  The  array  pattern,  in  each 
case,  for  n  =  200  can  be  regarded  as  the  best  simulated  opera¬ 
ting  pattern.  As  can  be  seen,  the  array  pattern  in  all  cases 
has  essentially  the  final  "best"  form  for  n  >_  8. 

All  the  results  in  Fig.  45  through  Fig.  49  were 
obtained  by  using  one  steering  vector  in  the  direction  of  the 
assumed  signal  direction  of  1°. 


6.2  Behavior  of  Output  Time  Axes  and  SNR  for  n  ■*  8 

Using  8  consequentive  samples  directly  ahead  of  the 
"current- value",  on  a  sliding  window  basis,  the  time-axis 
behavior  in  the  vicinity  of  the  "matched  point",  at  which 
acquisition  should  occur,  is  shown  in  Fig.  50  and  Fig.  51.  In 
Fig.  50,  the  time-axis  output  in  the  vicinity  of  8  consequetive 
peaks  corresponding  to  signal  epochs  is  portrayed  for  the  six 
steering  vectors,  with  one  interferer  of  N^  ■*  10^  at  67°,  a 
signal  of  S  =  1W  at  15°,  and  N  *=  o.lW.  The  ideal  SNR  ,  using 

A  © 

the  given  data  in  Eq.  7a,  was  found  to  be  SNR  ~  380.  The 

o 

simulation  run,  in  the  steering  direction  of  15°,  produced  an 
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SNRo  of  208,  calculated  by  dividing  the  peak  value  squared  by  the 
computed  variance  of  the  off-peak  values.  Thus,  we  were  able  to 
demonstrate  operation  within  2.6  dB  of  expected,  in  this  particu¬ 
lar  (set  of  8)  trials.  We  see,  then,  that,  for  this  particular 
simulation  an  adaptive  processing  gain  of  about  (30-2.6)  =  27.4  dB 
has  been  demonstrated;  if  the  24  dB  associated  with  the  255-long 
PN  processing  is  added,  a  total  of  about  51  dB  gain  has  been 
demonstrated. 

Fig.  51  shows  the  similar  results  if  now  two  500W 
interferers  at  67°  and  30°  are  present,  with  other  conditions  the 
same  as  Fig.  50.  Essentially,  this  demonstrates  operation  with 
multiple  interferers.  The  ideal  SNR  was  found  to  be  163;  this 
time,  during  the  particular  8  trials,  an  SNRq t 209  was  realized, 
which  exceeds  the  predicted  expected  value.  This,  of  course,  can 
happen  in  any  small-trial  case  when  the  SNRQ  of  concern  is  actually 
a  random  variable. 

To  emphasize,  in  viewing  Fig.  50  and  Fig.  51  it  must  be 
remembered  that  only  the  time  axis  in  the  vicinity  of  the  matched 
points  is  included,  to  give  a  representative  sample  of  non-peak 
behavior.  The  remaining  portions  were  exised. 
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Fig.  50.  Time  axis  output  for  8  signal  epochs  from  simulation  with  one  Interferer  of 
103W  at  67°  and  S  ■  1W  at  15°;  for  6  outputs  corresponding  to  six  steering 
vectors,  showing  behavior  only  in  vicinity  of  peak. 


output  for  8  signal  epochs  from  simulation,  with  two  Interferers  of  500W 
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7 .  Summary 

This  section  complements  the  Abstract  and  Summary 
(p.  i.)  by  reviewing  the  basic  issues  in  light  of  the  results 
described  above. 

A  particular  problem  in  space-time  processing  is 
the  enhancement  of  the  signal  radiated  from  a  "signal"  point 
source  at  unknown  distance  and  direction,  and  the  suppression 
of  one  or  more  signals  radiated  from  "interferer"  point  sources, 
together  with  the  supppresssion  of  directionless  "ambient"  ncise. 
The  use  of  multi-element  point  receivers  (an  array)  for  direc¬ 
tional  selectivity,  and  the  use  of  coded  signals  for  source 
selectivity  are  two  of  the  major  design  tools  in  creating  sys¬ 
tems  that  will  have  high  space-time  processing  gain. 

It  is  well  known  that,  for  mathmatical  "optimization", 
one  should  treat  both  time  and  spatial  processing  as  an  inte¬ 
gral  nondif ferentiated  processing.  However,  the  two  proceses 
are  separated  here  for  ease  of  implementation. 

The  spatial  processing  operates  to  reduce  the  system 
gain  to  each  and  every  point  source,  making  no  use  of  the  time 
character  of  such  sources;  of  course  the  time  character  is 
usually  unknown.  The  time  processing  uses  the  time  character 
of  the  "signal"  source  to  distinguish  it  from  the  other  souces. 
Something  has  to  be  done  to  prevent  the  space  processing  from 
reducing  the  gain  to  the  "signal"  source.  Thus  the  tradition¬ 
al  style  system,  with  the  spatial-processing  first,  is  forced 
to  violate  the  separation  of  time-and  space-processing  or  be 
doomed  to  defeat.  The  time-processing-first  system  proposed 
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and  discussed  here  attempts  to  overcome  this  problem  by  sort¬ 
ing  out  the  "signal"  source  first  (using  a  matched  filter) , 
and  passing  waveforms  on  to  the  space-processor  in  separated 
sets,  one  set  containing  nothing  from  the  "signal"  source, 
and  a  smaller  set  containing  waveforms  from  the  "signal" 
source,  and  a  smaller  set  containing  waveforms  from  all  sources 
but  with  the  signal  source  greatly  enhanced.  The  idea  is  that 
the  space  processor  "tunes  up"  on  the  waveform  set  that  does 
not  contain  the  signal  source,  and  then  applies  that  tuning  to 
the  set  of  waveforms  containing  the  signal  as  well,  thereby 
reducing  the  system  gain  to  the  non-signal  sources,  but  retain¬ 
ing  the  gain  to  the  signal  source.  However,  the  time-process¬ 
ing  can  only  do  the  separation  in  an  iterative  hypothesis-test¬ 
ing  fashion;  i.e.,  it  asks  the  subsequent  space-processor 
"what  would  the  output  be  if  this  set  were  signal  free  and 
this  other  set  were  to  have  the  signal  in  it?”  The  actual  op¬ 
eration  achieves  the  gain  of  a  joint  time  and  space  decision 
maker,  and  reflects  the  added  complexity  of  joint  time-space 
search. 

The  advantage  of  the  time-processing-first  configura¬ 
tion  is  that  it  is  an  open-loop.  It  is  a  well  known  folk  theorem 
that  open  loop  systems  have  a  lower  performance  threshold  than 
corresponding  closed  loop  systems;  that  is,  at  low  values  of 
whatever  limits  their  performance,  usually  some  SNR  or  power 
ratio,  the  open  loop  system  will  continue  to  perform  well  while 
the  closed  loop  system  fails. 

The  effort  here  has  demonstrated  the  expected  perform- 
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ance  of  the  PN  matched  filter  followed  by  an  SMI  spatial  pro¬ 
cessing.  Steady-state  performance  has  been  emphasized.  This 
initaial  evaluation  was  limited  to  studying  a  4-element  linear 
array  with  A/2  spacing,  and  consisted  first  of  writing  the  classi¬ 
cal  matrix  equations  for  the  performance  of  the  SMI  algorithm 
with  multiple  apriori  steering  vectors.  We  defined  the  array 
performance  measure  in  terms  of  the  output  signal-to-noise 
ratio,  SNRQ,  under  the  assumption  that  the  steering  vector  is 
pointed  in  the  direction  of  each  angle  considered.  The  input 
quantities  for  this  array  performance-measure  calculation  were 
the  matched-filter-output  signal  and  ambient  noise  powers,  and 
the  interferer  powers  and  directions.  These  quantities  were 
used  to  compute  a  deterministic  noise  covariance  matrix,  which 
was  then  inverted  and  multiplied  by  the  particular  steering 
vector  needed  for  each  angle. 

This  optimum  SNK0  performance  measure  versus  angle  for 
the  array  processing  was  programmed  using  BASIC  language,  for  use 
on  our  own  internal  minicomputer.  A  variety  of  interferer  cases 
were  run  and  studied  to  ascertain  optimum  performance. 

An  altered  version  of  this  program  permitted  studying 
the  theoretical  output  SNRQversus  angle  when  only  a  small  number 
of  steering  vectors  are  used.  This  permitted  studying  the  rela¬ 
tion  between  the  nulls  from  the  steering  vector  and  those  from  the 
interferers.  It  was  concluded  that  if  4  parallel  steering  vectors 
are  used,  along  with  a  decision  algorithm  that  responds  to  the 
highest  SNRQ,  one  can  provide  SNRQ  versus  angle  behavior  that  is 
nearly  omni-directional  except  for  the  angles  immediately  on  or 
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adjacent  to  the  interferer  angles. 

Another  large  computer  program  was  written,  using 
FORTRAN  language,  and  run  on  the  large  time-shared  computer,  that 
demonstrated  the  entire  operation  studied  here,  including  the 
matched  filter,  covariance  matrix  estimation  and  inversion,  and 
SMI  processing.  We  used  this  program,  for  a  few  interferer  cases, 
to  demonstrate  the  entire  operation  under  realistic  conditions 
where  the  covariance  matrix  must  be  estimated  from  arriving  matched- 
filter  output  samples. 

In  this  simulation  a  PN  matched-filter  of  length  255  was 
used  for  the  temporal  processing,  and  a  4-element  linear  array 
spaced  at  X/2  was  used.  Six  apriori  steering  vectors  were  used, 
placed  at  angles  of  +  15°  +  45°  and  +  75°  from  broadside.  The 
adaptive  array  output  for  each  of  these  was  computed  on  a  sliding 
window  basis.  The  outputs  are  compared,  with  a  threshold,  looking 
for  signal  peaks.  An  estimated  output  SNR  (based  on  a  limited 
number  of  noise  samples)  was  also  calculated  for  each  steering 
vector. 

The  desired  signal  consisted  of  2  or  more  periods  of 
a  255  bit  maximal  sequence.  This  signal  was  assumed  to  phase 
modulate  a  UHF  carrier.  The  source  was  assumed  to  be  located  at 
+15°,  which  corresponds  to  one  of  the  apriori  steering  vectors. 

The  phase  of  the  received  signal  at  each  antenna  element  was 
calculated.  The  computer  was  used  to  generate  the  complex  samples 
of  the  demodulated  signal  for  each  antenna  at  the  rare  of  4  complex 
samples  per  sequence  bit.  Similar  complex  samples  for  the  inter¬ 
ference  and  ambient  noise  were  added  to  the  signal  samples. 
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f./C. 
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The  matched  filter  outputs  for  filters  consisting  of  one 
period  of  the  sequence  were  calculated  for  the  in-phase  and  quadra¬ 
ture  component  at  each  antenna  element  on  a  sliding  window  basis. 

The  window  moves  in  increments  of  a  sequence  bit  to  avoid  missing 
a  signal  peak.  The  resulting  set  of  complex  numbers  are  the  signal 
inputs  to  the  adaptive  array  processing. 

One  or  two  interfering  sources  were  simulated  by  random 
number  generators.  They  were  assumed  to  transmit  carriers  phase 
and  amplitude  modulated  by  Gaussian  moise.  For  assumed  interference 
source  powers  and  source  positions,  the  received  amplitudes  and 
phases  were  calculated.  Corresponding  random  in-phase  and  quadrature 
demodulated  interference  samples  were  generated  in  the  computer  and 
added  to  the  corresponding  received  signal  samples.  For  each  an¬ 
tenna  element,  independent  in-phase  and  quadrative  samples  of 
Gaussian  noise,  representing  ambient  noise,  were  also  generated 
and  added  to  the  signal  samples. 

The  signal- free  complex  noise  samples  for  caluclating 
each  sample  covariance  matrix  were  obtained  from  the  outputs  of 
the  matched  filter  with  the  window  displaced  from  the  correct 
matched  signal  output.  From  1  to  200  samples  of  the  covariance 
matrix  were  calculated  and  averaged  to  investigate  the  effect  of 
the  number  of  samples  used  in  the  estimate  of  the  covariance  matrix. 

The  estimated  covariance  matrix  was  obtained  by  averag¬ 
ing  the  sample  matrices.  It  was  then  inverted  and  multiplied  by 
one  of  the  apriori  steering  vectors  to  find  one  set  of  combining 
weights.  The  matched  filter  signal  output  was  combined  using  these 
weights  as  one  possible  signal  output.  This  process  was  repeated 
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for  each  of  the  6  assumed  steering  vectors.  The  complete  process 
from  finding  a  new  matrix  to  calculating  6  more  possible  signal 
outputs  was  repeated  continously  on  a  sliding  window  basis. 

With  this  simulation  we  verified  the  fact  that  one  needs 
only  about  2k  =  8  conspijuetive  samples  to  estimate  the  covariance 
matrix.  The  resulting  output  demonstrated  processing  that  was 
within  3  dB  of  predicted,  since  interferers  of  103w  and  a  signal 
of  one  watt  was  used,  an  array  gain  of  30  dB  was  expected.  The 
simulation  resulted  in  an  array  gain  of  27  dB  which,  along  with 
the 24  dB  gain  of  the  matched  filter,  results  in  a  demonstrated 
57  dB  of  anti-interference  processing  gain. 

The  general  conclusion  of  the  above  predecessor  effort 
is  that  using  the  method  evolved  here,  including  the  4  parallel 
steering  vectors,  may  likely  be  the  best  way  to  combine  PN  waveform 
processing  and  adaptive  array  processing.  The  advantages  are: 

1)  the  combination  does  not  initially  (before  sync)  attempt  to  null 
the  signal,  but  offers  immediately  the  full  combined  processing  gain 

2)  the  response  to  either  signal  turn-on  or  interferer  variation  is 
near  instantaneous,  and  does  not  depend  on  loop  time-constants; 

3)  there  is  no  constraining  relation  between  interferer  size  and 
time -response,  as  there  is  for  a  closed-loop  system;  4)  the  de- 
correlation  of  wide-band  interferers,  at  array  input,  is  signifi¬ 
cantly  reduced;  5)  there  is  no  "great  race"  problem  between  modem 
synchronization  and  array  nulling,  as  there  is  when  a  PN  reference 
is  used  in  a  closed-loop;  and  6)  the  technique  here  should  be  more 
resistant  to  specific  and  determined  interferer  attacks  (such  as 
blinking) . 


It  is  recommended  that  the  potential  exhibited  by  these 
results  be  exploited  by  further  study.  In  particular,  it  is  recom 
mended  that  the  performance  evaluation  begun  here  be  generalized 
and  refined,  and  that  the  implementation  aspects  of  this  technique 


be  determined. 
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APPENDIX  A 

Review  of  Linear  Beamformer  Notation  and  Theory  Notation 

by  T.  G.  Birdsall 


X)  Let  *  denote  conjugate  and  T  represent  vector  or 
matrix  transpose 


2)  Let 

y  or  y(t)  be 

the  column  vector 

Of 

instantaneous 

inputs,  or  the 

complex  demodulation  of 

the 

inputs . 

* 

yT 

"  *yi'y2'  *'*  yk* 

(A.  1) 

3)  In  the  absence  of  signal,  i.e.,  when  the  input  is 
ambient  noise  plus  interferers,  let  M  be  the  autocovariance 
matrix  of  the  input.  We  assume  M  is  positive  definite 

M  =  E(yy*T)  *=  [E(yiyJT)]  (A.  2) 

M  is  hermetian,  i.e., 

M*T  »  M,  (M*-1)*7  »  M_1  (A. 3) 

4)  In  the  absence  of  interferers,  i.e.,  just  in 
ambient  noise,  we  assume 

M  -  al  (A. 4) 

5)  A  linear  beamformer  uses  weights  to  obtain  a 
single  number  for  each  "steering  direction" 

w  “  (W^,Wj,  ...  Wk)  (A. 5) 


•w: 
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For  a  given  set  of  weights  the  beamformer  output  Z  is 
2  -  W  *  y  =  W*Ty  (A. 6) 

Single-Look  SKEA+KGN  Theory 

The  theory  of  signal  detectability  provides  a  guide  in 
the  form  of  a  solution  if  the  signal  were  known  exactly 
except  for  amplitude  and  if  the  noise  were  wide-sense¬ 
stationary  gaussian  noise  with  zero  mean  and  known 
autocovariance  M.  It  states  that  if  V.  is  a  unit  vector 
proportional  to  the  input  if  it  were  caused  by  a  signal  from 
direction  q,  (in  the  absence  of  noise)  then  the  weights  to 
be  used  in  the  presence  of  noise  should  be  proportional  to 

W0  =  M^Vg  (A.  7) 

Performance  Evaluation: 

A  common  array  characterization  is  the  beam  pattern 
resulting  from  a  fixed  set  of  weights,  plotted  versus  the 
direction  of  arrival  of  a  unit  signal.  Let  0  be  the 
direction  used  in  setting  the  -weights,  and  <p  be  the  actual 
signal  direction.  * 

beam  pattern  =  WQ*TV^  =  (A.  8) 


The  beam  pattern  is  conjugate  symmetric  in  the  two  angles  0 
and  $,  and  gives  the  designer  or  evaluator  a  measure  of  beam 
width  and  sidelobe  levels. 


4L 

* 


58. 


More  direct  characterizations  of  the  beamformer 
performance  when  the  signal  direction  is  known  (or 
hypothesized)  would  be  the  beam  output  due  to  signal,  the 
beam  output  power  due  to  noise,  or  the  signal-to-noise 
ratio,  each  plotted  versus  the  direction  of  arrival  of  the 
desired  signal.  All  three  of  these  plots  are  identical  when 
the  weights 

W0  =  M_1Vq  (A- 7) 

are  used,  and  when  the  actual  direction  coincides  with  the 
steered  direction.  This  all-purpose  pattern  is  called  the 
PEP,  the  performance  evaluation  pattern,  PQ. 

pe  s  v0Tm“1v6  (A‘9> 

These  properties  will  be  quickly  proved. 


(a)  P q  «  beam  output  magnitude  for  unit  signal  in  steered 
direction 

(proof:  by  (A. 8)  P0  is  the  beam  output.  It 
is  real  and  positive  because  is  positive 
definite) 


(b)  P0  «=  mean-square  beam  noise  output. 

E ( |  Z |  ^  no  signal)  «=  E(ZZ*T) 


**r 

E(W  xyy  AW) 
*T  -1_ -1 


V6‘M 


MM  V, 


e 


*T 

W  MW 

*T  -1 
=  V6  M  A\ 


(proof : 
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(c)  Pfl  *=  beam  output  signal-to-noise  power  ratio 


(proof:  From  (a),  P.  is  the  beam  output 

magnitude  due  to  the  signal  in  the  steered 

direction,  so  P  2  is  the  beam  output  signal 

power;  from  (b) ,  P  is  the  noise  power;  hence 

6 

SNI>  -  *2 -  pe  >• 


Spatial  Whitening 

Although  theoreticians  love  the  weights  given  by  (A. 7) 
because  of  the  triple  property  of  the  P.E.P.  Pg  ,  it  is  often 
more  practical  to  adjust  the  weights  so  that  the  noise  on 
the  beam  output  is  independent  of  the  steered  direction. 
This  is  known  as  spatial  whitening 


For  spatial  whitening 


M^V, 


*T  -1 

VM  ve 


(A. 10) 


These  weights  are  just  as  optimum  as  those  given  by  (A. 7) 
since  the  theory  only  said  the  weights  must  be  proportional 
to  those  of  (A. 7) .  If  the  weights  of  (A. 10)  are  used  the  3 
performance  statistics  are 


Beam  output  magnitude  for  unit  signal  in 

steered  direction  **  /pT 

•  0 

Beam  output  mean-squared  noise  *  1 
Beam  output  SNR  «  Pg 


A. 11) 
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Example:  Linear  Uniform  Array,  Single  Interferer  plus 
Ambient,  Narrowband 


Consider  a  k-element  array  of  uniformly  spaced 
receivers,  and  a  unit  arrival  from  radians  off  endfire, 

o  i  a  <  i 


, , ,  nd  cos9v 

tt)  =  Yt  (t  — - —  ) 


(A.  12) 


where  d  is  the  element  spacing  and  c  is  the  speed  of 


propagation. 


Model  the  narrowband  signal  as  a  complex  phaser 


ya(t)  =  9(t)ei<2«ct+*!t>>  (A.  13) 


ZitfcndcosO,  ndcosO, 

yn+1(t,  .  9(t-  J«  e3(2’fct - 1  ei*<t~c— »  (s.14) 


The  usual  "narrowband  assumption" . is  that  the  array  delays, 
ndcosx/c,  can  be  neglected  in  the  amplitude  and  phase 
modulation  arguments,  so 

y„,(t>  -  g(t)eJ<2''£ct+»(t>'e-i2'fcnd  cos0/o  ,».15) 


For  brevity  we  define  the  "electrical  phase* 


0  *=  2irfcd  cos0/c 


(A. 16) 


A. 


J 


(A. 17) 


yn(t)  - 


y1(t)e 


When  d^c/2fc,  (A. 16)  may  be  inverted  to  find  the  physical  0 
from  the  electrical  5  . 

All  this  was  done  so  that  Vg ,  "a  unit  vector  in 
(electrical)  direction  0",  could  be  made  explicit 


(A.  18) 


The  entire  arrival  from  physical  direction  e  and  electrical 
direction  e  is  the  vector  valued  time  function 

y(t)  =  /K  VQ  y^t)  (A. 19) 

where  VQ  specifies  the  space  effects  and  y, (t)  the  time- 
effects.  The  spatial  autocovariance  is 

M  -  E(y  (t)y*T(t) )  «  VqV*t.  E  |yx Ct)  J 2  k  (A. 20) 

For  convenience  we  denote  the  matrix  V  V  as  m 

9  0  e 

«6  -  Ve* 

For  ambiant  noise  of  power  a,  plus  a  simple  interference 
with  power  b  from  direction  e ,  the  input  autocovariance  is 


V#"* 


s-5e 

-j26 


-j(k-l)0 
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M  =  al  +  bM 


6 


(A. 22) 


This  is  a  very  simple  autocovariance.  We  point  out  that 

lve!2  -  1  *m8  -  vev8TveveT  “  Vb’  “  Me  (A'23) 


Mq  does  not  have  an  inverse,  but  with  the  ambient  noise,  M 


(A. 24) 


does  have  an  inverse 

"'1  -  rl[r  -  S5E  Me] 

The  weights  for  steering  (signal)  direction  <p  are 

%  -  K'\  “  »_1K  -  m  v7%]  (A-25) 


Now 


V*TV  -  i  ky1  =  e3k(®"*)/2  sin (0-*)  )/sin(|(0-$) )  (A.26) 

.♦  kn =0  2  2 


This  is  the  well-known  beam  pattern  for  the  uniform  line 
array  Let  us  give  it  a  symbol 


Def 


B(x)  =  ejkx/2  sin(™)/k  sin(|) 


(A. 27) 


So 


Vqtv^  =  B(e-i) 


(A. 28) 
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w*  “  a'1[v*  "  i?b  B<e-*>ve]  lA-29) 


This 

has  the  intuitive 

sense 

of  pointing  toward  <J> 

and  away 

from 

e. 

The  performance 

evaluation  pattern  is 

P  =  V*TW 
*  *  4» 

_  “1 
-  a 

[x  -  5?E  B«e-*>vJTve] 

(A.  30) 

p* ' 

-  b 
a+b 

B(8-$)B($-0)] 

(A. 31) 

p*'a'1[1 

_ b_ 

a+b 

|  B(9-i)|2] 

(A.  32) 

Let  us  summarize  the  example,  using  a  non-unit  sig'va?  . 
Let  the  signal  have  power  s,  and  come  from  direction  )  let 
the  ambient  have  power  a,  and  the  single  interference  come 
from  direction  0  with  power  b.  The  output  signal  to  noise 
ratio  is 

SNR  =  (f>  "  s!e  I2] 

and  the  steering  weights  are  proportional  to 

»*  '  v*  -  iTE  Bl®-»)ve 


(A.  34) 
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APPENDIX  B 


Computer  Program  by  T.  G.  Birdsall 


•'edit  senate  p=j 

'jcomp  scri!aty-5-r2.j+2no  cnuto  p=j 

'sn d  cnatu 

DIM  X(l5),Y(l5),>J(lj),v(15)  'INPUT=X+JY,  0J7PU f=U+JV,  MATRICES 
DIM  56(4),  T6(4), 37(4), i/(4)  'I..RJT  i  JPEEHLid  V3C7 JR5 
PI=iPl/150  \  P2=2*iPI 

so  pr i n t\ph i ntnpr i nt\pr i ;i r  «c.mat5  <tgdi5aprilbi>" 

PRINT  "4  Ei.2ME.4i'  UNIFORM  Li  JS  A  Air  AY  (  .3  >«AVHL=NJTH  SPACING" 

PRiNT  " CM Afu  Cu.iPui'fi  BEAM  PATTERNS  WITH  J JI5E— ;iOri.«ED  STEERING" 
Co=.5  1 

FOR  £=J  fO  15  \  X(E)=£  \  Y(E)=3  \  NEXT  E 

PRINT  “ I NTcrtr ErtsisCE  AU70C0RRSLATI ON  MATRIX,  R,  IS  COMPLEX, 

PRINT  "HERmETI AN ,  i  T0PLIT2"  ) 

PRINT  M  THE  SUM  Or  V  7~*T  FROM  STRONG  DIRECTIONAL  SIGNALS,  AND11 
PRINT  «  A  UNIT  IDENTITY  DJE  TO  NON -DIRECTIONAL  AMBIENT  NOISE." 
PRINT  "THE  STEERING  VECTOR  Or*  RECEPTION  S  (SIGNAL  ONLY)  IS"  ) 
PRINT  "  ri  =  R- INVERSE-CONJUGATE  3" 

PRINT  "THE  PERFORMANCE  PATTERN  IS  N  DOT  5  FOR  ALL  DIRECTIONS  S" 
PRINT  "THE  WORMED  STEEdli^G  IS  N  /  SQROOTU  DOT  S)" 

PRINTXPRiMT  "AMBIENT  LEVEL  =  I  AT  EACH  SENSOR,  I.I.D." 

FOR  E=J  TO  15  STEP  5  V.  X(E)=I  \  NEXT  E 

PRINT  "PLEASE  DESCRIBE  THE  DIRECTIONAL  INTERFERENCE"  .« 

70  PRINT  "NUMBER  OF  SOURCES  ="J\INPUT  N7 
FOrt  1=1  TO  N7 

PRINT  "POWER,  DIRECTION  ANGLE  *  "(\INPUT  J,  A 
A1=P2*CO*COS(P1*A)  \  A2=A1/Pl 

FOR  R=1  TO  4  \  FOR  C=l  TO  4  \  E*4*R+C-5  \  A3=(C-R)*AI 
X ( E )  =X ( E )  + J*COS ( A3 )  \  Y(E)=Y(5)+J*SIN(A3)  \  NEXT  C  \  NEXT  R 
NEXT  I 

100  FOR  E-J  TO  15  \  U(E)=0  \  V(E)=3  \  NEXT  E  'SET  U+J7= IDENTITY  > 
FOR  E=0  TO  l o  STEP  5  \  U(E)=l  \  NEXT  E 

•'BIG  STEP  IS  To  INVERT  R,  THEN  CONJUGATE  RESULT 
'1ST  PHASE,  FORM  LOWER  TRIANGLE  ZEROS 
200  FOR  R=1  TO  4 

'NORMALIZE  SO  DIAG  ELEMENT  =1  ) 

E=3*R-5  \  N=H  \  A=X(E)  \  B=Y(E)  \  G0SU3  600 
'CREATE  ZEROS  BELOW  THAT  1 

M=R  \  FOR  N=R+1  TO  4  \  E=4*N+M-5  \A=X(E)  \3=Y(E)  \G03U3  5f*J  \  N.  IT 
NEXT  H 

'2ND  PHASE,  FOriA  UPPER  TRIANGLE  ZEROS 

FOR  C=4  TO  2  STEP  -1  \  FOR  R=1  TO  C-l  \  E=4*R+C-5  ) 

N=R  \  M=C  \  A=X ( E )  \  B=Y(E)  \  GOSUB  503 

NEXT  H  \  NEXT  C  \ 

'NO*i  CONJUGATE  RESULT  j 

FOR  E=0  TO  15  \  V (E)=-7(E)  \  NEXT  E 

3oJ  PRINTXPrtINr  "STEERING  ANGLE  START,  STOP,  STEP  ="J\INPJT  Z1,Z2,Z. 
PRINT  "CALC  ANGLE  START,  STOP,  STEP  »«*\IMPJT  .1  ,W2,..3\PRINT 
FOR  A=Z I  TO  Z2  31HP  Z3 
PrilA'i  "JTEErIN J  A.<GLE  =  "/»}A 

print  "jNrjrmed  weight  vector*" 

PRUT  "REAL","I..iAj AG", "PHASE  (DEG)" 

GjGO^  4.g  '37+JT/  IS  .(EIG.;.'S,  RB  AT  Ai3  IS  PERrOR.iANCE 
FOR  1=1  TJ  4  \  a ?-6 / ( I )  S  Y?=r/(I)  \  GOSUB  3300  \  AV=3oOwAV  ) 

PRINT  M8.a5*XV,YP,RV,  N6.D2JA?  \  NEXT  I  j 


I 
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PRINT  “PSArod-ikACZ  ...EASUHE  =  “TS.GoJR  ii  “  AT  "  ;  A6.  }2 ;  A3 ;  '•  DEJ;*  \ 

P.U.i’T  "TOR.iED  .iSIJ.iT  VcCTJ.R  USED  Iif  BEAM  F.'ifTE X.i  IS" 

HR  I  .<i  *' .-i E AL “  *  "  I  -iAJ  •'  y  “MAG"  ( •VriASE" 

.-ijs  jJrt(  .io)  1 

Pod  I = I  TJ  4  \  S7 ( I )=37( I  )/Rd  \  T7( I )=T7( I )/R3 
X9=S  /  <  I  >  \  Yv=T7  ( I )  \  JOSJJ  3666  \  A9=36D*A9 

PRINT  Z8.0:>*X9,Y9»R9,  *6.J2tA9  \  TEXT  I  ) 

HfU.IT  “ - (BEAM  PD<.ER  PATTErT) - " 

PR  ITT  "aTOLIE"  ,  "POWER  G*IT" 

FOR  3=<i  I  TJ  «2  STEP  R3  \  B1*P2*CJ*C05CP1*3>  ' 

FOR  1=1  TJ  4  \  Sot  I >=C03 ( (I- 1 )«a I )  \  To( I )=SIT  ( ( I-l)*31 )  \  NEXT  I 
GJSUd  6  JO  'CALC  THE  DOT  PkJDJCT 

PRINT  Xi By  T8.03;R3*R3  ) 

TEXT  3  \  Ha I NT  \  NEXT  A 

39J  PRITT  "l=>..iJRE  ATGLE3,  2=>CHAHGE  INTERFERENCE  *  MS\INPUT  0 

IF  0=1  THEN  30J  \  IF  0=2  TiiET  5.)  \  GOTO  390  .) 

'5J34JJ  FIT  jo  THE  nEIJiiTS  A  OT-BEAM  DOT  PRODUCT,  THE  PERFORMANCE  EVALUAT.' 

4JJ  AI=P2*CJ  *Cg5(PI*A)  \  A2=Al/HI  'A  I  IS  RAD,  A2  13  DEJ,  (ELEC  PH  A'.  ) 
FOR  1  =  1  TO  4  \  Sod)  =C0S  ( ( I  - 1 )  -V  A I  )  \  To  (  I )  =S  I T  ( ( I  - 1 )  *A  1 )  \  HE  XT  i 
JOS  Jo  106  \  003  JB  S  36  'FIND  STEERING  VECTOR  AiJD  DOT  PRODUCT 
RETJrtT  ,) 

'SUjoJJ  R Jr  A  =  RG.i  if  -  (A+Jo)  *  RO.i  .4  -  N=/=M 

oJJ  FOR  C3= I  TJ  4  \  £l=4*T+Co-o  \  E2=4*M+C5-5  ") 

X(EI  )  =X  (  £  l ) -A*X ( E2 ) ( E2 ) 

Y  ( E  l ) =Y ( E I )  -A*Y ( E2 ) -G*X ( E2 ) 

U(El)=U(cl)  -  A*  J  ( £2 )  -s-osv  V  ( E2 )  ) 

V ( E I )  =  V ( £  I ) - A*7 ( H2  > - J * J ( £2 ) 

NEXT  Co  \  RETURN 

) 

'SUooOJ  DI/IJE  ROW  .J  BY  A+Jb 

oJJ  Z6=A*A+is+3  \  Ao=A/XJ  \  B6=-3/Z6  'Ao-J-JBo  =  1/  A+JB 

FOR  C6=l  TO  4  \  E=4:fii+Co-5  ) 

'OLD  VALUE  =  SJ-5-J  To,  TER  VALUE  =  S7*J  T7 

56=X(E)  \  To=Y(E)  \  30SU3  660  \  X(E)=S7  \  Y(E)=T7 

S5=J(E)  \  To=V(E)  \  G05U3  666  \  U(E)=S7  \  V(S)=T7  ) 

Te*T  Co  \  RETURT 

004  S7s.\6*36—36*76  \  T/=A6*76+tS5*36  \  RETURN 

) 

'303/ JO  MATRIX  U+J  V  4  VECTOR  S6+J  To  =  VECTOR  S7+J  T7 

y  VECTOR  So+J  To  IS  UNCHARGED 

7JJ  FOR  R=l  TO  *r  \  Rl=4*R-o  \  S/(R)=J  \  T7(R)=D  )  '■ 

FOR  C=!  TO  4  \  E=R  1  i-C 

3  7 ( R)  =  57 ( R) * J ( E ) *3 j(C )-V ( 2 ) *To( C) 

i /(R)  =  T7(:t)  +  V(E).vSo(C)+U(E)-.-To(C)  ) 

TEXT  C  \  id  XT  R  \  RET'urN 

'SJ^o  ..'j  X9--J  Yv  =  VECTJR  (So-J  id)  jJT  (S'/+J  T  / )  ) 

'  RE  i  j  J  jy  AS  R  j  Ai  AS  DEGR  EES 

oJJ  :<y=;  \  i'ys=v  \  FJ.t  i  =  l  TO  4 

X/  =  Xy+Jo(  I)...,  /(!)  +  .  j(I);ri  /(I) 

Ty='yyJo(I)..-i /d)-To(I)  .'3/(1)  \  TEXT  I 
uOJJJ  3 ^  >4.'  \  \  Ao=3oJ-'/.y  \  ..ETUitT 

t 

'Cj.«i  i.i'JE  .ii  T.i  id/ycRD  3Jo3  .  .  I.«  R2H 


) 
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